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ABSTRACT 


After a general introduction, in which some features of the density distribution and 
velocity distribution are discussed in connection with the possibility that the galaxy 
may have a spiral-like structure, the author attempts in the third section to elucidate 
the relation between density distribution and distribution of peculiar velocities. It is 
shown that if the large-scale unevenness in density found ina recent study is accepted, it 
is natural to expect a deviation of the stream vertex from the longitude of the galactic 
center of similar order as the deviation actually observed. 

The distribution of light in the nearly spherical nebula NGC 4494 and in the strongly 
flattened elliptical nebula NGC 3115 has been investigated with the aid of plates taken 
by Dr. Oosterhoff with the 60-inch reflector on Mount Wilson. The results, as shown in 
Table 1 and Figures 3 and 4, can be fairly well represented by Hubble’s formula 
IT = Io/(1 + r/a)?, though some small but probably real deviations are indicated. The 
isophotes in NGC 3115, of which a few typical ones are shown in Figure 5, deviate 
markedly from ellipses and show remarkably sharp points at the equator. The ratio of 
the axes varies with the distance from the center; it has a minimum for a semimajor 
axis of 1’ (see Fig. 6). The space distribution of light in NGC 3115 has been calculated; 
a section with isophotal surfaces is shown in Figure 7, while the run of space density in 
the equatorial plane may be found in the last column of Table 1 and in the dotted curve 
in Figure 4. The light density near the center is of the same order as that in the center 
of the globular cluster M 3. 

In the last section these data are considered in relation to Humason’s measures of 
the rotation of NGC 3115. It is found that beyond about 10” from the center in the 
equatorial plane the mass density must be about constant. The actual density depends 
on the unknown flattening of the attracting mass, but it should be at least of the order of 
140 solar masses per cubic parsec (10~?° gm/cm3), or 2000 times that near the sun. In 
this region of constant mass density the light-density diminishes with a factor of at 


t This paper was presented at the symposium on galactic and extragalactic struc- 
ture, held in connection with the dedication of the McDonald Observatory on May 5-8, 
1939. 
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least 10; the distribution of mass in the system appears to bear almost no relation 
to that of light. The total mass of the system cannot yet be estimated; it should be at 
least 5.10'° solar masses but may be much larger. Expressing both mass and light in the 
sun as unit, one finds the ratio of mass to light in the outer parts of NGC 3115 to be 
about 250; this contrasts sharply with the value of 1.8 found in the neighborhood of the 
sun and suggests an extreme relative density of very faint dwarfs or interstellar ma- 
terial. The analysis further indicates that the peculiar velocities of the material (stellar 
or other) responsible for the light of the nebula must be high, between 100 and 200 
km/sec for the parts within @ = 20’’, between 175 and 300 km/sec near a = 45”, the 
actual values depending again primarily upon the flattening assumed for the attracting 
mass. The distribution of velocities perpendicular to the equatorial plane indicates a 
remarkably strong deviation from a normal Gaussian law (cf. Table 5). The results 
of this section should be regarded as quite tentative, the main uncertainty being due to 
the possible existence of important absorption and diffusion of light within the nebula. 
I. THE GALACTIC SYSTEM 
STRUCTURE 

About the beginning of this century the outstanding problem for 
the explorer of the universe was that of the size and general shape of 
the galactic system. The Plan of Selected Areas initiated by Kap- 
teyn in 1906 was proposed principally with a view to this problem. 
The Plan has been successful in that it has helped to obtain a fairly 
reliable picture of the general features of the density distribution 
within a sphere of, roughly, 1500 parsecs radius around the sun. On 
account of the strong absorption of light near the galactic plane it 
has not yet been possible to make much use of it in outlining the 
dimensions of the system in the plane of the Milky Way, which, 
however, could be derived from other data, notably those on the 
distribution of globular clusters, and recently also of cluster-type 
variables, gathered by Shapley and others. 

Though it can hardly be said that our knowledge concerning the 
general shape of the system is anything like final, yet some of the in- 
terest is shifting gradually to another phase, namely, that of detect- 
ing an eventual spiral structure. As has been pointed out by several 
authors, it can hardly be doubted that in all observable character- 
istics (such as extreme flatness, occurrence of many highly luminous 
stars, existence of a well-defined center) the galactic system resem- 
bles most strongly that class of extragalactic systems which invari- 
ably show a more or less well-defined spiral structure. 

The view that the galaxy might have a spiral structure has been 
expressed almost since the first discoveries of spiral structure in 
nebulae. The oldest explicit reference to this seems to be in a paper 
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by Stephen Alexander.? Both he and Proctor,’ seventeen years 
later, tried to find support for their spiral theory and to construct to 
some extent the galactic spiral from the appearance of the Milky 
Way as inferred mainly from the observations of the Herschels. In 
1900, and more completely again in 1913, Easton,‘ apparently quite 
independently and unaware of these earlier suggestions, made a very 
careful study of the Milky Way as shown by modern visual and 
photographic observations, especially with a view to delineating the 
course of the possible spiral arms of the galactic system; he placed 
the center of the spiral in the direction of Cygnus. 

Later developments have made it clear that this representation 
can hardly resemble the real structure of the galaxy, the main cause 
of its failure and of the failure of all attempts to find spiral-like 
structure being the strong and uneven absorption near the galactic 
plane; for there is little doubt that many of the features of the appar- 
ent Milky Way structure are determined rather by the distribution 
of absorbing material than by that of the stars. 

With the data at present available the difficulty (which, as indi- 
cated, is essentially that of separating the effects of absorption and 
of real structure) can be overcome more or less satisfactorily only in 
regions at some distance from the Milky Way, where the total 
amount of absorption may be estimated from the numbers of extra- 
galactic nebulae (and, in some cases, also from color excesses). The 
procedure involves two assumptions, which it has been possible to 
verify only to a certain extent: (1) the assumption that the distribu- 
tion of faint nebulae presents no large-scale irregularities and (2) 
the supposition that the absorbing material is mainly confined to re- 
gions within a few hundred parsecs from the galactic plane. The 
available evidence on these two points is still unsatisfactory. 

2 “On the Origin of the Forms and the Present Condition of Some of the Clusters of 
Stars and Several of the Nebulae,” A.J., 2, ror, 1852. 

I am indebted to Dr. Rosseland for this reference as well as for that to Proctor’s 
paper. Dr. Rosseland has also drawn my attention to the possibility that Lord Rosse 


himself may have had the idea that the galaxy had a spiral structure, but so far no 
direct published reference to this by Lord Rosse has been found. 

3 “A New Theory of the Milky Way,” .N., 30, 50, 1869. 

4“A New Theory of the Milky Way,” Ap. J., 12, 136, 1900, and “A Photographic 
Chart of the Milky Way and the Spiral Theory of the Galactic System,” ibid., 37, 105, 
IQ13. 
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The results of a study of the Selected Areas along these lines, if 
accepted at their face value, show that the sun is situated in a rela- 
tively sparse region, surrounded in practically all the longitudes 
covered by the observations by large regions where the star density 
within 500 parsecs of the galactic plane is higher than that near the 
sun. These results are certainly very tentative, and it is probable 
that more extensive and more adequate subsequent studies may in 
several regions yield considerably different densities. But the general 
trend of the density distribution would seem to be fairly well estab- 
lished. Naturally, then, the thought presents itself that the density 
variations found bear some relation to the spiral structures observed 
in extragalactic systems; but the indications are still vague and en- 
tirely qualitative. 

In a flat system like the galaxy, the main results of such density 
analyses as the one described can be most easily represented by com- 
paring the counted numbers of stars in a given region to the number 
which would be found in a standard system in which the equidensity 
layers are planes parallel to the galactic plane and are so arranged 
that they conform to the density distribution perpendicular to the 
galactic plane as found in the vicinity of the sun.’ The true density 
distribution in space can later be easily derived from the residuals 
found. I believe that in the present preliminary stage of this type of 
investigation the procedure of starting with this simple comparison 
has considerable advantage, because it does not involve the use of 
any other data than those directly observed, and at the same time it 
brings out the essentials of the density phenomena studied. 

It has been found by several astronomers that the arms in spiral 
nebulae have a color corresponding to a fairly early spectral class.° A 
large-scale survey of faint B- and A-type stars in our galaxy, accom- 
panied by a systematic and precise investigation of their colors, 
might accordingly furnish a good opportunity for a further search for 
such a structure in the galactic system. Radial velocities and studies 
of differential galactic rotation would furnish a necessary check on 


5 See Vashakidse, Bull. Abastumani Obs., No. 1, p. 87, 1937, and No. 2, p. 109, 1938; 
Oort, Ann. d’ap., 1, 71, 1938, and B.A.N., No. 308, 1938. 

6 Seares, Proc. Nat. Acad., 2, 553, 1916, and Carpenter, Pub. A.S.P., 43, 204, 1931; 
compare also an unpublished article by Seyfert on relative colors of arms and nuclei. 
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the average total absorptions. Beside color excesses and galactic- 
rotation effects, proper motions of faint stars of known spectral type 
may become an important help in evaluating the absorption. 


MOTIONS 

It may be in place here to mention a few points concerning stellar 
motions, which would seem to have a bearing on the problem of the 
internal structure of the galactic system. 

So far, the most important information derived from proper mo- 
tions of faint stars has been concerned with the velocity distribution, 
as it has been shown that this distribution remains practically identi- 
cal over a large region in the galactic plane, extending to about 1500 
parsecs in the northern Milky Way. Not only do the motions of 
stars over this entire region follow the same ellipsoidal law, but they 
also appear to show the same lengths for the axes of the velocity 
ellipsoid and thus to indicate that the “temperature’’ of the stellar 
system does not vary perceptibly throughout the region investigated. 

The extremely valuable study of radial velocities of faint 6 Ceph- 
ei variables recently published by Joy’ adds important evidence to 
that already in existence on the relative motions of different parts 
of the galactic system. All available data concerning these motions 
can be adequately interpreted as indicating that the motions of the 
various parts are closely perpendicular to the direction of the mate- 
rial center of the system, as found from globular clusters (J = 325° 
+ 2° p.e.) and independently also from the distribution of stars at 
large distances from the galactic plane® (J = 324° + 3°). From a 
comprehensive study of the “rotation-effects” in the radial velocities 
Plaskett and Pearce? found / = 325°4 + 1°1 (p.e.); Joy’ obtained 
5°3 + 1°3 (p.e.), while Merrill and Sanford’s recent study of 
velocities from interstellar lines'® yielded 329°. The asymmetry of 
high velocities and our systematic motion with respect to globular 
clusters and near-by extragalactic systems also point in a direction 
roughly go° away from the center. Considering these data together, 
there appears to be pretty strong evidence that the motions in the 

7 Ap. J., 89, 356, 1939; Mt. W. Contr., No. 607. 


§ B.A.N., No. 308, p. 263, 1938. 9 Victoria Pub., 5, 281, 1936. 
10 Ap. J., 87, 118, 1938; Mt. W. Conitr., No. 585. 
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galactic system, at least in our surroundings, are circular around the 
center. If systematic motions along spiral arms exist, the arms must 
be very nearly circular. The data for differential rotation now extend 
over a disk of nearly 2500 parsecs radius around the sun, except for 
longitudes between 193° and 271°, where only a few observations are 
available. Apart from a (probably) local effect in the bright, rela- 
tively near-by B-type stars, the data show no evidence of a syste- 
matic expansion or contraction; Joy, for instance, finds a quite negli- 
gible value of —3.8 km/sec for the K-term in his distant Cepheids. 

In the bright Cygnus cloud Joy has observed two faint Cepheids, 
GL and V343 Cygni, which, judging from unpublished color meas- 
ures recently obtained by Oosterhoff, appear to be situated at dis- 
tances of about 6000 and 10000 parsecs and are yet not more than 
350 and 600 parsecs from the galactic plane; MZ Cygni, at —10° 
latitude, likewise seems to have suffered but moderate absorption 
and to be about 5000 parsecs distant. If Joy’s velocities for these 
three stars may be considered as representative for the regions in 
which they are situated, they show that also in this very distant 
part of the galactic system the motions can be satisfactorily repre- 
sented by a simple rotation of this system. Computing the distances 
from the radial velocities, on the assumption that the residual veloci- 
ties are wholly due to rotation of the galaxy, I find for GL Cygni 
6700, for V343 Cygni g1oo, and for MZ Cygni 5500 parsecs, in very 
good agreement with the values just quoted, which were derived 
from the apparent magnitudes. 

The available data seem as yet to be unsuited for a significant 
determination of the distance to the center and still less of the 
arrangement of mass within the galactic system. Berman," and, 
more recently, Camm” have published evidence which they inter- 
preted as indicating that the smallest planetary nebulae are situated 
at distances comparable to that of the galactic center, and some even 
considerably farther. They conclude that the velocities of these 
nebulae give, therefore, information about the dynamics of the sys- 
tem in regions near the center. The evidence, however, is not con- 
vincing.*3 

tt Lick Obs. Bull., 18, 57, 1937. 


12 M.N., 99, 71, 1938. 13 Cf. Oort, M.N., 99, 376, 1939. 
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The large-scale irregularities indicated in the density distribution 
do not seem to modify the gravitational field to any observable ex- 
tent, at least there is no sign of this in the differential motions. This 
need not cause surprise, for it can be inferred from the available data 
that as much as 85 per cent of the total gravitational force in our 
neighborhood may be considered as due to a central mass and as 
varying approximately as the inverse square of the distance from the 
center, while the rest can be adequately ascribed to a homogeneous 
mass extending beyond the sun. 


THF DEVIATION OF THE VERTEX 


Though the systematic motions of stars with small peculiar mo- 
tions may be almost uninfluenced by local unevenness in the density, 
circumstances are quite different with regard to the peculiar veloc- 
ities. It can easily be seen that a close relationship must exist be- 
tween the distribution of peculiar velocities and the distribution of 
density in the surroundings of the sun. To simplify, let us consider 
stars moving exactly in the galactic plane. Size and shape of a stellar 
orbit may then be characterized by two quantities, for which we may 
take the mean distance from the center, @,, and the velocity II, with 
which the star will pass the circle with radius @, (these parameters 
are evidently equivalent to the parameters 7, and J2, the energy in- 
tegrals and the momentum integrals, used so extensively by Lind- 
blad). As long as there is no instability and as long as I], is fairly 
small compared to the rotational velocity of the system, the motion 
of a star can be described approximately as following an elliptical 
epicycle, the center of which moves with uniform velocity over a 
circle of radius @, around the center of the system." Suppose, for 
simplicity, that the frequency distribution of II, is the same over the 
whole region considered and also that the phases of the motions of 
the stars along the epicycles are distributed uniformly. The distribu- 
tion of the centers of the epicycles will then evidently be directly re- 
lated to the density distribution, while it also determines the veloc- 
ity distribution at any given point. 

It is well known that in a system having rotational symmetry the 

'4 This description was first used by Lindblad in 1927; cf. Arkiv. mat., astr. 0. fysik, 
20A, No. 17, and Upsala Medd., No. 26. 
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resulting distribution of velocities should be symmetrical with respect 
to an axis pointing to the center of the system, and this axis should 
be the major axis of the velocity ellipsoid. 

The observations, especially those of bright stars, show, however, 
that the longitude of the major axis of the velocity ellipsoid deviates 
between 15° and 20° from the direction of the center. For the faint 
stars a deviation in the same direction is indicated, although there 
seems to be some evidence that it would be less pronounced for these 
stars. 

In the representation of the system just given the most direct in- 
terpretation of the vertex deviation is that it indicates a lack of rota- 
tional symmetry in the density distribution. The assumption of an 
asymmetry in the density distribution of the right order to explain 
the deviation of the vertex is quite in line with the unevenness in the 
density discussed in the first section. Inversely, it might be said 
that the existence of the vertex deviation indicates large-scale den- 
sity deviations. The relation between the two can best be illustrated 
by considering a simple hypothetical case. 

Consider the motion of a star whose velocity is small compared to 
the rotational velocity of the galactic system. Let w be the distance 
from the center of the system and w, the mean distance. If 0 is the 
transverse velocity and 0, the same velocity at @, then the equa- 
tions of motion can, in first approximation, be written as 


° ‘ ‘s , ® —w)O0K; , 8 
Ow = O.a., 6 = K+ o — &,)0l —+—., 
0a w 
Restricting ourselves again to first-order terms in (@ — @,)/@o, we 
can put 
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so that the second equation may be reduced to 


& = — 0(@ — w), 


where 


02 Vo be 
n= 2 (24. a 3) — 


@ Wo Ow 


A and B denoting the usual constants of differential galactic rota- 
tion. We obtain, accordingly, 


®—@ =asinn, 


while the difference between the star’s actual angular velocity, w, 
and the angular velocity, w., of the center of the epicycle described 
by the star, is evidently given by 


0 0, 24: Os... 
O-W = -—--—-= —=--snn. 
Ww Wo Wo Wo 


In a system of co-ordinates rotating with an angular velocity w., 
the orbit of the star is, therefore, given by 


x=oO—@ =asinn, 
2d 

y = — (A — B) cos nt, 

: n 


y being the transverse co-ordinate reckoned from the center of the 
epicycle in the direction of the rotation. 


Assuming A = +0.018 km/sec-parsec, B = —o.013 km/sec - 
yarsec,'® and, accordingly, 7 = +0.040 km/sec-parsec, we obtain 
I , $1) 4 Pp 


an orbit as shown in Figure 1. Let us now consider a point S situated 
in a longitude differing by ¢ from that of the galactic center. We are 


'S M.N., 99, 374, 1939: 








interested in the “peculiar velocity” at S or, more precisely, in the 
difference between the original star’s velocity II, 6 in the point S and 
the velocity corresponding to circular motions at S. Denoting this 
circular velocity by 0.,s, we have 
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0.8 = 0. + 


and 


so that 


while 


II = an cos nt = 2By. 


If a is the angle between the 
peculiar velocity at S and the 
direction of the galactic center, 


we thus have 


(8 — @,5) 


tga = —~— 


Suppose now, for an illustra- 
tion, that, as seen from the sun, 
the density of centers of epicycles between 248° and 271° galactic longi- 
tude and from 1400 to 1900 parsecs distance were twice that in the 
corresponding longitudes at the other side of the galactic center (be- 
tween 19 and 42°). At an average longitude of 260° (g = 115°) the 
distance limits considered would correspond to x = 560 and 730 
an cos nt = 2Bx tg ¢ = —0.056x, 
those stars near the sun whose epicycle centers lie in the region con- 
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sidered would have peculiar velocities in the II-direction between 33 
and 44 km/sec. With the aid of the foregoing equation, we find that 
the angles a corresponding to the limits of g are 193° and 216°, so 
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Fic. 2.—Velocity ellipsoids. The numbers attached to the ellipses indicate frequen- 
cies relative to that in the center. 


that the double star density in the region considered entails a dou- 
bling of the frequency of peculiar velocities in the neighborhood of 
the sun between, roughly, 35 and 50 km/sec in position angles be- 
tween 193° and 216° as reckoned from the direction of the galactic 
center. The corresponding changes in the velocity distribution are 
roughly pictured in Figure 2. Ordinates are peculiar velocities in the 
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II-direction, abscissas those in the 9-direction. The full-drawn el- 
lipses, with their major axes pointing to the galactic center, represent 
lines of equal frequency in a symmetrical stellar system; the accom- 
panying numbers indicate the relative frequencies expressed in terms 
of that at the center of the figure as unity. The local increase of star 
density considered causes the lines of equal frequency to bulge out, as 
shown roughly by the dotted curves, the bulges having been smoothed 
so as to avoid sudden transitions. The bulge extends over the range 
of curves with frequencies between about 0.40 and 0.16. A similar 
bulge for the inner curves may be obtained by an increase of the 
star density between 700 and 1400 parsecs in the same longitude in- 
terval by 50 per cent. Though we cannot make a very certain esti- 
mate of the vertex deviation which would be caused by the bulges 
considered, it is probable that it would come out at least of the same 
order as that observed; at the same time, it need not entail a very 
noticeable effect on the systematic motion of the stars considered. 

The foregoing example does not pretend to resemble the actual 
conditions in the galactic system. It is evident that a relative excess 
of density in a sector opposite to the one considered would produce 
the same deviation in the vertex. It is also possible. and even proba- 
ble, that conditions are complicated by the fact that any unevenness 
in the density distribution is likely to be accompanied by an uneven- 
ness in the distribution of the phases of the motions. 

Available density data are still insufficient to indicate whether 
density deviations in the required sense exist; all that can be said, I 
believe, is that the data do not contradict the interpretation sug- 
gested. The only purpose of the present discussion was to illustrate 
the relation between velocity distribution and density distribution 
and to indicate the order of the deviations in the density required to 
explain the deviation observed in the velocity distribution. It may 
be well to point out, in this connection, that rather different inter- 
pretations are possible, such as, for instance, that proposed by Lind- 
blad."° He suggests that the orbits of the stars responsible for the 
deviation of the vertex extend to a region where circular motions are 
unstable. The orbits of these stars would then be of a different char- 
acter from those considered, starting in a spiral-like manner in the 


16 Stockholms Obs. Ann., 12, No. 4, 1936; Osten Bergstrand Festskrift, 1938. 
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unstable region; it is easily conceivable how, on this model, a slight 
difference could occur between the velocities of outgoing and of in- 
coming stars, which would give rise to a deviation of the vertex. 

Lindblad has pointed out that the more or less sudden appearance 
of the strong asymmetry in the velocity distribution which is found 
to occur at peculiar velocities of about 65 km/sec can also find ex- 
planation in his model. The most direct interpretation of the sudden 
appearance of the asymmetry would, however, seem to be that it is 
caused by the way in which the density in the galactic system varies 
in its outer regions. If considered in the same way that we consid- 
ered the vertex deviation, it would mean that, after remaining rough- 
ly constant up to about 1600 parsecs in the direction of the anti- 
center, the density would rapidly fall off beyond this distance. Such 
a variation could certainly be reconciled with the available data. 


II. EXTRAGALACTIC SYSTEMS 

Accurate knowledge of the density distribution in extragalactic 
systems is evidently one of the prerequisites for a study of the dy- 
namics of these systems. Apart from a few color measurements, as 
reported above, little or no quantitative work has yet been done on 
spiral nebulae proper. For the elliptical nebulae quite interesting re- 
sults have been obtained by Hubble, and more recently studies on 
the same subject have been published by Redman and Shirley. A 
number of these nebulae have been studied in some detail in Leiden, 
with the aid of plates taken by Dr. Oosterhoff at the Mount Wilson 
Observatory and by the author at the Perkins Observatory. The 
measurements and reductions were carried out in collaboration with 
Mr. Kleibrink. In the following pages a report will be given on the 
results for two of the nebulae, NGC 4494 and NGC 3115. The re- 
sults given are provisional in several ways, and some further meas- 
urements and reductions are contemplated. 

The most interesting outcome of Hubble’s original investigation 
of the elliptical nebulae’? was that, except for differences in flatten- 
ing, all these nebulae appeared to be built on the same model. In 
each of the 15 objects studied the surface intensity could be well rep- 
resented by the formula J = J,/(1 + r/a)*; and this is not all: it ap- 


17 Ap. J., 71, 231, 1930; Mt. W. Contr., No. 398. 


ee ee 


ee se 


Se Sete on eR IES 


Prt epee 


Hee a oe 








286 J. H. OORT 


peared also that the “central” intensity, or rather, the parameter J,, 
was nearly identical for the various nebulae; the same may hold ap- 
proximately for the linear diameters. No satisfactory reason has yet 
been found why these nebulae should all be built on this particular 
model, but it is evidently a result of the greatest interest.** 


PHOTOMETRY OF NGC 4404 

Figure 3 shows detailed results of the intensity measures on NGC 
4494 as studied at Leiden. Four plates taken by Dr. Oosterhoff with 
the Mount Wilson 60-inch reflector were used, with altogether six 
exposures ranging from 1 to 60 minutes. The nebula had been classi- 
fied as Eo by Hubble, but the measures indicate a ratio of 0.83 for 
the shortest and the longest axes. For the plot in Figure 3 various 
position angles were simply averaged (for most exposures tracings 
were made through 12 different position angles). The ordinates in 
Figure 3 show logarithms of surface intensity /, this being expressed 
in stars of magnitude 25.0 per square second of arc as unit.'? The 
intensity scale was determined by means of sensitometer marks; on 
two of the plates it was also determined from measurements of sur- 
face brightness in the rings of far interfocal images of a Selected 
Area exposure superimposed on the field of the nebula. These Se- 
lected Area exposures were also used for obtaining the zero points of 
surface brightness. The two independent determinations of the in- 
tensity scale agree very well. A partial third test is furnished by the 
intercomparison of the various exposures. Zero-point corrections 
have been applied to the different exposures in order to fit the curves 
together. It will be seen from the diagram that, with these correc- 
tions applied, the various exposures agree quite satisfactorily. Sig- 
nificant differences are found only for points derived from the faint 

8 Tt should be mentioned that Redman (M..N., 96, 588, 1936; see also ibid., 98, 613, 
1938) has recently found somewhat deviating results for some of the nebulae measured 
by Hubble. On the other hand, Miss Patterson, working at Harvard, has in an unpub- 
lished study obtained data which generally confirmed Hubble’s earlier measures. The 
results for NGC 4494 and 3115 can be represented quite well by Hubble’s formula. 

19 The surface brightness expressed in magnitudes per square second of arc (//) is 
related to J by the formula H = 2.5 (10 — log J). The relation between the surface 
intensity 7’, expressed as the number of stars of 10.0 photographic magnitude per square 
degree (as has also been used in recent work) and that used in the present note, /, is 


I’ = 13.01. 
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and ill-defined part of the toe of the reduction curve. Observations 
corresponding to these uncertain parts have been indicated by re- 
duced size of the characters in the figure. Abscissas in the graph are 
logarithms of the distance to the center in seconds of arc. The sur- 
face brightnesses plotted have not yet been corrected for the effect of 
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Fic. 3.—Distribution of light in NGC 4494. Ordinates are logarithms of surface 
intensity, expressed in stars of 25™o per square second of arc. Abscissas are logarithms 
of the distance from the center in seconds. The dotted curve represents Hubble’s 
formula log J = log /, — 2log (1 + r/a), with log J, = 3.41 and a = 3"0 (1” is rough- 
ly 12 parsecs). 


atmospheric and instrumental diffusion of the light; this must have 
had a pronounced effect upon the curve to the left of log r = 0.5, but 
it is believed to have been of relatively little importance for the rest 
of the curve.”° 

20 Redman and Shirley have recently shown (M.N., 98, 614, 1938) that on their 


plates a measurable fraction of the light is dispersed to distances of several minutes from 
the source of light. Qualitatively, this has been confirmed on the Mount Wilson plates 
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As will be seen from the dotted curve in Figure 3 as well as from 
the numbers in Table 1, the main features of the density distribution 
in this nebula (which has not been measured by Hubble) can be 
quite well represented by Hubble’s formula, with a = 3%0 and 


TABLE 1 


DISTRIBUTION OF SURFACE BRIGHTNESS IN NGC 4494 AND 3115 


























NGC 4404 NGC 3115 
r LOG r Major Axis | Minor Axis | , 
Major 
log Iobs*| O-—CtT Axis 
| 
llog Tobs*| O-Cf |log Jobs*} O—Ct | log ot 
| | 
B-Olsascdiise — © 3.28 | —o 13| 3.61 | +o0.26] 3.61 | +0.06 +3.26 
see rae 0.0 3.20 | + .04| 3.48 | + 21| 3.46 | + 17] +2.91 
“aN NE eee ae 0.4 2.92 | + .04/ 3.18] + 02] 3.07 | + .07| +2.19 
4.0 0.6 2.65 | — .03) 3.03] — .03| 2.83] + .o4] +1.83 
Ot eee 0.8 2.38 | — .05| 2.89 | — .04] 2.56] + .03} 41.51 
Oe ree 1.0 2.1% | — .03] 2.73 | — .03] 2-94 | oo} +1.19 
ORS Nona tes 1.2 1.83 | + .o1| 2.55 | + .o1] 1.85 | — .07| +0.87 
re eee 1.4 OsS0 ) th 204) 3.30) [Og F487 — .o9g} +0.48 
IS oils ae 1.6 1.16] + .06] 2.04] + .07/ 1.16 | — .04] +0.01 
OE Se ge 1.8 0.78 | + .o5|} 1.71 | + .07/ 0.85 | + .03| —o.50 
BOI cs ce cas 2.0 0.34 0.00) 1.29] + .oI] 0.51 | +0.08) —1.06 
158.0. . ae Be ; | 0.86 | — .06) | -1.65 
to: ee | Bi fa cnscaeficees ss] Oigg | 0.99 —2.27 
| 





* The symbols J and ¢ are expressed in stars of photographic magnitude 25.0 per square, or cubic, 
second of arc as unit. 

t The columns O—C contain the differences between the observed distributions and those computed 
from the formula log J = log J. — 2 log (1 +7/a). The constants used are NGC 4494, log Jo = 3.41; 
a = 30; NGC 3115 major axis, log Jo = 3.35, @ = 102; NGC 3115 minor axis, log 3.55, a 2784 

t The last column shows the distribution of light in space along the equatorial plane. 


log J, = 3.41. There are, however, some minor deviations which are 
probably real; between 4” and 8” the computed values are about 0.06, 
or o™15, brighter than the observed intensities, while between 25” 


used in this article by tracings through the stars BD-+26°2359 (84, Ao) situated only 
6’ from NGC 4494, and BD—6° 3066 (10™2 pg, K2), 19’ from NGC 3115. The magni- 
tude of the effect may be illustrated by stating that 3 per cent of the light is dispersed to 
distances between 3’ and 2’ from the star. The percentage is found to be practically the 
same on all four long-exposure plates measured; it appears to be independent of seeing 
conditions and might well be due to dispersion in the plate. The amount is about four 
times less than that found by Redman and Shirley at the same angular distances. The 
effect upon the outer contours is small, even for the very condensed nebula NGC 4494. 
In order to correct for it, the value of log 7 for r = 100” should be diminished by 0.04, 


that for r = 63” by about 0.02. 
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and 50”’ they are about the same amount fainter. It should be men- 
tioned that, on account of the uncertainty as to the true distribution 
of light near the center, the parameters a and J, cannot be sep- 
arated with much accuracy. A value of a = 4"0, for instance (with 
log J, = 3.20), would fit the well-observed part of the intensity 
distribution quite as well as the values adopted. 

The distance of NGC 4494 may be estimated from its radial veloci- 
ty, which was kindly communicated to me by Humason, in advance 
of publication. As corrections to the center of the galactic system 
are negligible in this case, the distance can at once be obtained by 
the aid of the velocity-distance relation. Using Hubble’s value of 
530 km/sec per million parsecs™ for the coefficient of expansion and 
Humason’s velocity of +1350 km/sec, we find 2.5X10° parsecs. 
The nebula is only 14° away from the nucleus of the Virgo cluster; if 
it is a member of the nuclear part of this cluster, its distance should 
be 2.2 10° parsecs.””7 Assuming, provisionally, a distance of 2.4X 10, 
parsecs, 1’’ at the distance of NGC 4494 is equivalent to 12 parsecs, 
and we find a = 36 parsecs. 


PHOTOMETRY OF NGC 3115 

Of NGC 3115, three 60-inch Mount Wilson plates taken likewise 
by Dr. Oosterhoff were measured, containing altogether five ex- 
posures of 60, 503, 15, 4, and 1 minutes, respectively. Details of 
measurements and reductions will be more fully published elsewhere 
in combination with the results for the other nebulae studied; we 
shall discuss only the main results of this study. 

Several circumstances appear to make this elliptical nebula espe- 
cially interesting for thorough investigation. The object is the bright- 
est of the nebulae classified as E7 by Hubble. Its flattening, corre- 
sponding to a ratio of about 3 to 1 between the two axes, is the 
strongest observed among the nebulae classified as elliptical; this 
indicates that it is an extreme case, which, in the current thoughts on 
spiral formation, may mark the point of transition to an unstable 
system; it also indicates that the nebula must be seen practically 
edge on, and this knowledge enables us to derive the space distribu- 
tion from the distribution of light as seen projected on the sky. Fur- 


2 The Realm of the Nebulae, p. 170. 22 Tbid., p. 163. 


Ce oe ee ee 


Gee ees ee eee 








290 J. H. OORT 


thermore, the rotation has been measured by Humason, the meas- 
urements extending to about 45” from the center. 

Figure 4 and Table 1 show the results of the measures of the dis- 
tribution of light in the nebula as projected on the sky, along the 
major and minor axes, respectively. The various exposures have 
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Fic. 4.—Distribution of light in NGC 3115. The full-drawn curves show the varia- 
tion of surface intensity J along the major and minor axes. The 1-minute exposure is 
represented by crosses, the 4-minute by dots, the 15-minute by triangles, the 50}-minute 
by circles, and the 60-minute by squares. The dashed curve shows the variation of the 
space density of light, y, along the equatorial plane. Abscissas are logarithms of the 
distance from the center in seconds, ordinates are logarithms of the light-intensity ex- 
pressed in stars of 250 per square second for J (numbers on the left) and per cubic 
second for g (numbers on the right). 


again been fitted together by applying zero-point corrections. In 
this case a provisional correction has been applied in the region near 
the nucleus for the dispersion of light caused by atmosphere, instru- 
ment, and plate. This dispersion, as determined from registrograms 
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through a focal star image 80’’ from the center of the nebula on the 
one-minute exposure, is +0.04 mm on the plate, or +170. The cor- 
rections applied to log J range from about +0.18 in the center to 
—o.08 at 3” distance, becoming negligible beyond about 4”. No 
claim is made for the numerical accuracy of these corrections; the 
most that can be said is that they indicate the order of the uncer- 
tainty arising from the light dispersion. It would seem that the shape 
of the isophotal contours (as shown in Fig. 5) and the results for the 
concentration of the light near the major axis (except in the central 
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Fic. 5.—Isophotes in NGC 3115. The isophotes correspond to log J 2.4, 2.0, and 1.6, 
I being expressed in stars of 25™o per square second. The co-ordinates are given in 
seconds (1” is roughly 7 parsecs). The 1-minute exposure is represented by asterisks, 
the 4-minute by crosses, the 15-minute by squares, the 50}-minute by dots, and the 
60-minute by circles. 


nucleus) can have been but little disturbed by the small observa- 
tional dispersion. 

The distributions in Table 1 can again be fairly well represented 
by Hubble’s simple formula, though there are some deviations which 
appear to be real. The differences from the values computed by this 
formula, using Hubble’s results for the parameters a (viz., 1072 for 
the major axis and 284 for the minor axis), are as shown in Table 1; 
the values of log 7, used for fitting the formula to the present ob- 
servations are 3.35 for the major axis and 3.55 for the minor. 

Isophotes for log J = 2.4, 2.0, and 1.6 are shown in Figure 5. 
They were obtained from a considerable number of tracings in di- 
rections parallel with and perpendicular to the principal axis. Spe- 
cial care was taken to obtain tracings near the equatorial plane and 
accurately parallel to it, in order to get as reliable data as possible on 
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the run of the isophotes near the major axis; for this purpose each 
exposure was traced along the equatorial plane as well as at 0.07 and 
0.20 mm (or 19 and 5"4) distance on either side of it. The exact 
situation of the major axis had been ascertained by micrometer 
measures and had been indicated on the plates by dots; in this way 
the error of centering in the microphotometer has been reduced to 
something of the order of 071, while the probable error of the orienta- 
tion will not have exceeded o°1. The various exposures have been 
indicated by separate signs; the corresponding isophotes were re- 
duced to exactly the same size by applying factors derived by equat- 
ing the average of the minor and major axes of the isophotes con- 
sidered to the standard values of the axes for the surface brightness 
concerned as read from Figure 4. 

The most remarkable features of the isophotes are the remarkable 
sharpness with which they run into the major axis and the long, 
straight portions which lead down to these sharp edges—features 
that have also previously been remarked upon by Hubble;?s his data, 
however, show a little less detail in the parts near the axis. There is 
an indication that for the inner isophote, with semimajor axis a = 
21’’, and especially for isophotes still nearer to the center, the sharp- 
ness of the edge disappears, while at the same time the ratio c/a of 
minor to major axis increases. The ratio reaches a minimum of 0.29 
near ad = 60”; it increases again for the outer isophotes. The varia- 
tion of c/a with log a between a = 4” and a = 300” is shown in 
Figure 6, in which the ratios for various exposures have been plotted 
separately. Both the innermost isophotes (say for a below 10’’) and 
the outermost isophotes (beyond a = 100”) are somewhat uncer- 
tain on account of the instrumental and atmospheric dispersion of 
light. The effect upon the outer isophotes has been determined on 
the long-exposure plates 460 and 467 by measures on the near-by 
tenth-magnitude star BD — 6° 3066 (cf. also n. 20, p. 287); the 
curve gives the corrected values, while the crosses have been plotted 
without correction. In order to strengthen the part of the curve for 
large values of a, a number of registrograms were made of a fine 100- 
minute exposure taken by Ritchey on December 25, 1911, which 
was kindly put at my disposal during a stay at the Mount Wilson 


23 Ap. J., 71, 264, 1930; Mt. W. Contr., No. 308. 
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Observatory in the summer of 1939. The results, again uncorrected 
for dispersion, have been plotted as asterisks in Figure 6. Hubble’s 
earlier determinations of the axial ratio*4 for the brighter parts of the 
nebula have also been indicated in the figure (as squares); they are 
slightly higher than the values derived here. No attempt has yet 
been made to correct the inner isophotes for the effect of light-dis- 


persion. 
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Fic. 6.—Axial ratio of isophotes in NGC 3115. Dots represent results from 1-minute 
exposure, circles from 4-minute, and crosses from 503- and 60-minute exposures. The 
asterisks have been derived from a 100-minute exposure by Ritchey. The squares have 
been copied from results of short 100-inch exposures previously published by Hubble. 
The curve has been approximately corrected for the effect of light-dispersion over 


great distances (see text). 


SPACE DISTRIBUTION OF THE LIGHT IN NGC 3115 

In order to obtain at least some rough insight into the true struc- 
ture and possible dynamics of NGC 3115, the admittedly provisional 
isophotes discussed in the foregoing section were used to determine 
the space density yg. Let z denote the distance from the equatorial 
plane of the system; let ¢,(@) define the space density of light at dis- 
tances @ and z from the minor and major axes, respectively, ex- 
pressed as the number of stars of 25.0 photographic magnitude per 
‘cubic second of arc” at the distance of the nebula. Let J,(x) be the 
surface brightness at a distance x from the minor, and z from the 


24 Ap.J., 71, 264, 1930; Mt. W. Contr., No. 398. 
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major axis. If it is assumed that the system considered has rotational 
symmetry around an axis perpendicular to the line of sight and that 
there is no absorption of light within the system, it can be shown 
that in any plane of constant z 











‘i dI,(x) 
g(a) = -< , dx dx .?5 
: x? _# @? 


aw 


In most of the actual integration, which was carried out by Mr. 
Pels, it proved easiest to determine first, graphically, the derivative 
of log J with respect to log x and then to compute d//dx from these 
quantities. The extrapolation over the unknown outer regions is of 
little importance. If, in rough accordance with the observations, it 
is assumed that for large values of x (say, beyond X) J becomes of 
the form cx~? and if this formula is assumed to hold to infinity, the 
part of the integral expression for g beyond X is easily seen to be 


c @ @? ._ W 
—s- 5 ~~ mail — = => ae ee. 
Tw ( X J a” = 


If X is taken to correspond to the outer isophotes measured, the ex- 
trapolated part is usually only a small fraction of the total. The cal- 
culations were carried out for 8 values of s and were roughly checked 
by the inverse calculation of J from ¢. 

Figure 7 shows a section through the isophotal surfaces as provi- 
sionally found from these calculations. The curves resemble, in gen- 
eral, the isophotes in the projected image, but the ratio between 
minor and major axis comes out somewhat larger; between a = 30” 
and a = 100”, c/a is found to be approximately 0.40, as against 0.30 
in the projected image. However, it should be remembered that the 
true flattening may well be somewhat stronger, as the line of sight 
may be slightly inclined to the true equatorial plane. The values of 
log yg corresponding to the various surfaces are indicated in the 
figure. In the center we find, approximately, log g = 3.26. Taking 

2s The formula was first given by Von Zeipel (Ann. de l’obs. de Paris, Mém., 25, 
p. F30, 1908); see also Smart, Stellar Dynamics, pp. 297-302; it is easily seen to apply 
to any surface of revolution the axis of which is perpendicular to the line of sight. 
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the photographic light of the sun per cubic parsec as unit of density 
of light?® and estimating the distance of the nebula as 13 million 
parsecs?’ (therefore, 1’’ = 7.3 parsecs and m — M = 25.9), we find 
that a star of 25™o per cubic second corresponds to 0.74 units; if the 
true distance is x times larger, the corresponding light per cubic 
parsec will be « times smaller. The central intensity corresponds to 
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Fic. 7.—Section of isophotal surfaces in NGC 3115. The numbers in the figure are 
logarithms of the light-density ¢ corresponding to the respective curves, ¢ being ex- 
pressed in stars of 25™o per cubic second. For the three inner curves log ¢ is, respec- 
tively, 2.33, 2.03, and 1.73; in the center log ¢ is estimated to be 3.27. Co-ordinates are 
in seconds (1” is roughly 7 parsecs). 


1400 units. This is 30,000 times the total density of light in the part 
of the galactic system in which the sun is situated and is of the same 
order as that in the center of Messier 3.”° 


* The photographic absolute magnitude of the sun was assumed to be +5.26, in 
accordance with a recent discussion by Kuiper (A p. J., 88, 429, 1938). 

27 This distance is still very uncertain. The only criteria from which it can be 
estimated are the radial velocity, the apparent magnitude, and the size. The radial 
velocity, corrected to the center of the galaxy, is +450 km/sec, corresponding to a 
distance of 0.8 million parsecs; however, the smallness of the velocity makes this result 
doubtful. From the apparent magnitude, 9.8 according to Shapley and Miss Ames, 
combined with Hubble’s value of —15.2 for the average absolute magnitude of this 
class of nebulae (The Realm of the Nebulae, p. 176), a distance of 1.0 million parsecs is 
derived. The apparent dimensions, on the other hand, indicate a much larger distance: 
the average of dmax and @min for NGC 3115 is 6"5, while the average linear value of a for 
the 7 elliptical nebulae measured by Hubble in the core of the Virgo cluster amounts to 
66 parsecs; if we suppose that the average a for NGC 3115 is the same, we find 2.1 
million parsecs for the distance. The foregoing value is thus a compromise; it should be 
kept in mind that the true distance may well differ from it by a factor as large as 2 in 
either direction. 


28 Hertzsprung, A.N., No. 4952, 1918. 
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The run of log ¢ with @ in the equatorial plane of the nebula has 
been plotted as a dotted curve in Figure 4. The curve for log ¢ runs 
considerably steeper than that for log 7; from ® = 80” on, ¢ can be 
expressed as w~3". The smallest light density recorded is approxi- 
mately g = 0.01, or 0.007 units as previously defined, or one-sixth of 
the light-density in the region surrounding the sun. 


DYNAMICS OF NGC 3115 
It was hoped, originally, that the distribution of the light in extra 
galactic nebulae would bear resemblance to the distribution of mass 
and that investigations of the light-distribution might give some 























TABLE 2 

‘ @ ® —Kz Oc Oo 9 log ¢ 
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co ee oe 234 28.9 2.70 44.1 316 2.36 
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insight into the question of how far the most flattened elliptical 
nebulae and the regular cores of early-type spirals approach un- 
stable configurations which could give rise to the formation of spiral 
arms. Humason’s results on the rotation of NGC 3115 seem to indi- 
cate, however, that there is not much prospect in this direction; not, 
at least, for this nebula. 

Table 2 shows the values of the potential ® and the acceleration 
K, for points in the equatorial plane of the nebula, computed on the 
assumption that the distribution of density is the same as that of the 
light. The factor by which the true distance is larger than the dis- 
tance of 1.5 million parsecs provisionally assumed is x; f is the un- 
known ratio of mass density to light-density, the former being ex- 
pressed in the mass of the sun per cubic parsec, the latter in the 
photographic light of the sun per cubic parsec as a unit (in the direct 
vicinity of the sun f is about 2). The potential was assumed to be 
zero at infinity. 
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The values of ® and Kz were obtained by numerical integration 
with the aid of the equidensity surfaces given in Figure 7, and with 
rough estimates about the more outlying surfaces. The surface cor- 
responding to log g = — 2.17 (@ = 233’) was arbitrarily taken as the 
boundary of the nebula, the density being supposed to become zero 
outside this surface. The neglect of the material outside this surface 
is not likely to have had an appreciable influence on K,, but the 
values given for the potential will be too low. Accepting the approxi- 
mate values given, the velocity of escape for © = 470 is found to be 
I 19 V «f km/sec; with « = 1 and f of the order of 100, it may rough- 
ly be estimated to be of the order of 1200 km/sec. 

It is evident from the run of the values of ® in Table 2 that, if the 
mass distribution in a system like NGC 3115 is assumed to be like 
that of the light, such a system would be very stable. The condition 
for the appearance of dynamical instability is that at some place the 
force varies more rapidly than @~;, or the potential more rapidly 
than *. With the numbers given in Table 2, the steepest variation 
of log ® with respect to log occurs in the outer parts of the nebula, 
but even at @ = 120” we have only ® ~ o@~°-**, that is, an even slow- 
er variation than in the case of a Newtonian field, so that the system 
would be very far from a state in which it could begin to develop 
spiral structure. 

However, the spectrographic measures of rotation indicate that 
the distribution of mass is probably very different from that of the 
light; in such a case the stability cannot, of course, be judged from 
the distribution of the light. Humason’s rotational velocities are 
given under 9, in the sixth column of Table 2.2? The measures repre- 
sent the average radial velocity of the various parts of the nebula 
seen projected on top of each other; it is easily seen, however, that 
in a case like the present, where the rotational velocity increases 
linearly with the distance to the center, the observed radial velocity 
will be equal to the true rotational velocity at the corresponding 
point. Humason’s measures extend to 45” from the center; on ac- 
count of the small scale the measured velocity for @ = 871 is likely 
to have a considerable percentage of uncertainty. 

It is interesting to compare the run of the rotational velocities with 


29 Report of the Director of Mt. W. Obs., 1936-37, p. 31. 
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the entirely different run of the “circular” velocities (velocities cor- 
responding to circular orbits) computed on the assumption that the 
mass distribution is the same as that of the light; these are given in 
the preceding column under 9, and appear to remain practically 
constant over the range in » considered. As 0, will, in general, be 
lower than 9,, we infer from the data for @ = 4575 that, if x = 1 
(that is, if the distance is 13 million parsecs), f must be at least of the 
order of 150. Furthermore, if f were taken constant, the circular 
velocity at 871 and 1671 would come out many times larger than the 
rotational velocity. On the other hand, if we should assume that the 
rotational velocities are of the same order as the circular velocities, 
we should find that f varies with » and that the distribution of the 
density differs radically from that of the light. 

I believe it can be shown that the first alternative is very improb- 
able. Let us assume for a moment that it were true; suppose, there- 
fore, that f be constant and, for instance, equal to 200.°° 

The general relation between 9, and 9, in systems which are in 
a dynamically steady state may be written" 


ag loge log IP 6 — 0.) | 
2 = 82 = Jf? ae g 108 ¢ =— g 108 ‘. ae (9 _— 85)’ , (1) 


dloga <dloga q Tr 


Here II and 9 denote the velocities in the direction of @ and in that of 
the rotation; bars denote average values, thus 8 = @,. In the case 
of NGC 3115, where the angular velocity of rotation is approximate- 
ly constant, the term {1 — (0 — 9,)?/II’} is probably small; in the 
accompanying tentative calculations it has been neglected. 

The seemingly perfect regularity and symmetry of the system 
suggests an approximately steady state. Jeans* has shown that in 
the case of strict dynamical equilibrium the velocity distribution in 
the II-direction must be identical with that in the Z-direction (paral- 
lel to the axis of rotation), so that 





a i. (2) 


3° The total mass of the system would then be 7 X 10'° solar masses. 

3 The equation can, for instance, be derived directly from equation (10) given by 
Jeans in M.N., 82, 122, 1922. 

32 See, for instance, M.N., 76, 81, 1915. 
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Although it is by no means clear how such an equality of the two dis- 
tributions would be established**—for during an interval of 10'° or 
10" years the transfer of energy between individual stars in NGC 
3115 must have been negligible—I shall nevertheless use the above 
equation as a tentative working hypothesis for the following calcula- 
tions. It should be noted that, as far as an estimate of the line width 
can be made, it indicates, according to information from Humason, 
that the dispersion of the II-velocities cannot be much larger than 
about 100 km/sec—an amount, therefore, which is of the same order 
as that computed for the velocities in the Z-direction (see Table 5). 

If the gravitational force is known, the distribution of velocities 
in the Z-direction can be derived from the density distribution. In 
the computations I have assumed that the distribution of velocities 
could be represented as the sum of a number of exponential distribu- 
tions, as follows: 


$2) = Dae 8, (3) 


in which }°@; = 1. We have then the following equation,*4 in 
which K, represents the gravitational acceleration 


g:(@) _ roe. i (4) 
Po(@) 


The mean square velocity for a point near the equatorial plane is 
given by 
‘a 0; 
} a ee 
al: (5) 
The computation of this mean square velocity might also have 
been made with the aid of the more general formula 


= 3 if ¢:K,dz, (6) 


A 
ye. 


‘I 


which, however, would have necessitated extrapolations to infinity. 


33 Tam indebted to Dr. Woltjer for a discussion on some of the difficulties involved in 
introducing this equation. Cf. B.A.N., No. 238, p. 253. 
34 Cf. ibid. 
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In the hypothetical case considered above, where the distribution 
of mass is the same as that of light and the ratio of mass to light is 
200, K, can be directly computed from the known light-distribution. 
Somewhat approximate* calculations were made for two regions of 
the nebula, namely, at @ = 1671 and 32” 2, respectively. The results 
are shown in Table 3 (units cm/sec? and cm?/sec’, respectively; x was 
assumed to be 1). 

From these data it was found that the velocity distribution for 
@ = 16"1 could be well represented as the sum of two Gaussian dis- 
tributions with mean square velocities of 160 and 335 km/sec and 
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with 6, = 0.52 and 9, = 0.48, so that the total mean square veloc- 
ity VZ = 259 km/sec. For & = 322 we obtain also a sum of two 
Gaussian distributions, corresponding to mean square velocities of 
95 and 280 km/sec with 6, = 0.60 and 6, = 0.40; the total mean 
square velocity is 192 km/sec. 

We can now insert these values of Z? in formulae (2) and (1) and 
compute 9, from the latter. As has been explained, the last term in 
the parentheses was neglected. From the values of Z? just found, 
d log II?/d log & was estimated to be about —0.86; 0 log ¢/4 log a, 
which has been computed from the data on the light-distribution, is 
given in the last column of Table 2. In this way we find 8, = 448 and 
468 km/sec for @ = 16%1 and 32" 2, respectively. These values differ 
very considerably from the values 676 and 620 km/sec computed 

3s The approximation consisted in replacing the real equidensity surfaces by sphe- 
roids, which were so chosen that their volumes as well as axial ratios equaled those of the 
corresponding equidensity surfaces. 
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from the light-distribution. Such high values for the circular veloc- 
ity could have been obtained only by using very much larger values 
(between 400 and 500 km/sec) for the mean square peculiar veloc- 
ities, values which would have been in contradiction with the dis- 
tribution of light in the z-direction. Probably such high values as 
400 or 500 km/sec for the mean peculiar velocities could likewise be 
dismissed on the ground that the widths of the spectral lines would 
then have had to be greater than those observed. It may be con- 
cluded that the distribution of mass in the system must be consider- 
ably different from that of the light. A study of the light-distribu- 
tion does not, therefore, enable us to draw conclusions regarding 
such problems as stability or the total mass of the system. 

The question now arises: How must the distribution of mass den- 
sity be in order that both the observed distribution of light and the 
measured rotational velocities as well as the conditions for dynami- 
cal equilibrium can be satisfied? The available data are not sufficient 
to permit a complete answer. We can start from attracting masses 
of various degrees of flattening. Detailed calculations have been 
made only for a model in which the equidensity surfaces in the at- 
tracting body are ellipsoids of axial ratio c/a = 2/3. The high value 
2/3 was selected in order to be safely above the limit 0.552, where the 
instability of a homogeneous mass sets in. The computations were 
made for three points, at @ = 8%1, 1671, and 45"5, respectively; the 
gravitating mass was supposed to consist of an ellipsoid shell I with 
semimajor axis d = 10”, an ellipsoidal shell II extending from a = 
10” toa = 20”, and another shell III from a = 20” to a = 70”, out- 
side of which the density was assumed to be zero. The densities in 
the inner ellipsoid and in the two shells were now determined in such 
a way that for the three points mentioned the equations (1), (2), 
and (4) were satisfied. The results are in Tables 4 and 5. 

The mass densities and light-densities are expressed in terms of 
the sun per cubic parsec. It will be noted that throughout the nebula 
the mass density is extremely high. Beyond a = 10” a value of 140 
solar masses per cubic parsec is found, roughly equivalent to 10~”° 
gm/cm:, or nearly 2000 times the density near the sun. The total 
mass of the system would become 5 X 10°° solar masses. As the extent 
of the outer ellipsoid is entirely arbitrary, the real mass may differ 
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considerably from this figure; the system may easily be as massive as 
the galactic system, the mass of which is estimated to be between 
10% and 2X10". 

The most remarkable feature is that the mass density appears to 
remain practically constant outside a 


TABLE 4* 


DISTRIBUTION OF MASS AND LIGHT IN NGC a115 
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tensity falls off steeply. The strongly condensed luminous system 
appears imbedded in a large and more or less homogeneous mass of 
great density. This is a feature which, to some extent, might be ex- 
pected at once from a superficial consideration of Humason’s re- 
sult that the nebula rotates like a solid body; for, if the peculiar 
velocities are moderate, so that, except near the center, the measured 
rotational velocity resembles the circular velocity, Humason’s re- 
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sult means that the attracting mass must approximately be a homo- 
geneous ellipsoid enveloping the measured region. 

It may be noted that the mean square velocities, as shown in the 
fourth line from the bottom of Table 5, are still high, though they 
are considerably smaller than in the model considered previously. A 
notable feature is the pronounced non-Gaussian character of the dis- 
tribution of the velocities Z, which suggests that we are looking at a 
mixture of very different types of stars, or of stars and interstellar 
material. Though the particular velocity distributions found depend 
to a considerable extent on the assumptions concerning the mass dis- 
tribution, we have seen that the non-Gaussian character remains, 
even if the mass distribution is supposed to be the same as the light- 
distribution. : 

A very rough calculation along the same lines as that described 
but assuming axial ratios of } instead of } for the attracting ellip- 
soids indicates that in this case also the densities in the outer ellip- 
soids would come out nearly equal; the absolute densities as well as 
the values of f, however, would become more than twice as high as 
those found previously, while the peculiar velocities would also in- 
crease by a factor of nearly 2. The mean square velocity at @ = 
45"5 would become of the order of 300 km/sec, a value so large that 
it could probably be verified from observations. 

In the outer parts of the nebula the ratio f of mass density to light- 
density is found to be very high; and this conclusion holds for what- 
ever dynamical model we consider. The spectrum of the nebula 
shows the characteristics of G-type dwarfs. Since f cannot be much 
larger than 1 for such stars, they can account for roughly only $ per 
cent of the mass; the remainder must mainly consist either of ex- 
tremely faint dwarfs having an average ratio of mass to light of 
about 200 to 1 or else of interstellar gas or dust. 

Let us first consider these latter two possibilities. If the mass con- 
sisted mainly of gas, the optical thickness for rays coming from the 
nucleus would be of the order of tooo times the average optical 
thickness traversed by the light of a star at 1000 parsecs distance in 
the galactic plane. The presence of layers of gas of such thickness 
might be verified from suitable spectra. 

Many objects of the apparently more or less related type of Sa 
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spirals show strong evidence of absorption effects. Though no direct 
trace of such effects can be found in NGC arts, the possibility can- 
not be excluded that there would be considerable amounts of solid 
interstellar material distributed evenly through the nebula. At the 
same time, however, it must be observed that if this matter is at all 
like the absorbing matter in the galactic system, it can at most ac- 
count for a negligible fraction of the total density in NGC 3115, 
which is at least 3000 times that of the interstellar material in the 
galaxy. For, if the great density in the nebula were due to such ab- 
sorbing matter, the transparency would be so low that we could not 
penetrate more than a very few parsecs into the system.*° 

The high value of f found stands in sharp contrast to what is 
found in the galactic system. In the vicinity of the sun f is about 2; 
it may rise to, at most, 7 at large distances from the galactic plane.*7 
The ratio for the bulk of the galactic system can hardly be very 
much larger. A minimum value for the total amount of light con- 
tained in the part within the sun’s distance from the center may be 
estimated by assuming that the light-density within an ellipsoid 
having a semimajor axis equal to the distance of the sun from the 
center (adopted as 8000 parsecs) and an axial ratio c/a = 1/20, is 
everywhere the same as that observed near the sun, i.e., 0.048 times 
the photographic light of the sun per cubic parsec. The total amount 
of light thus found is 0.51X10"° units. From the rotation of the 
galaxy the total mass within the sun’s distance from the center may 
be estimated as 9.1 X10'°, so that we should find f = 18. The true 
value of f must be expected to be considerably smaller, as the light 
will probably be concentrated toward the center, so that the above 
estimate would give too small a value for the total light. We may 
conclude that the ratio between mass and light for the galactic sys- 


36 Tt is conceivable that an absorption might exist which is sufficiently powerful to 
reduce the light of the stars in the system by a considerable factor and that, thus, the 
real ratios f for the undimmed light would be much smaller than those inferred from the 
present discussion. Such absorption would, however, also affect the central regions and, 
if called in to explain the high f-values in the outer region, would give impossibly low 
values for the central parts. 

37 Cf. B.A.N., No. 238, p. 285, 1932. It should be noted that the unit of light in that 
article was a star of +6.0 photographic absolute magnitude; the light densities given 
should therefore be multiplied by a factor of 0.506 to be reduced to the present units. 
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tem as a whole is probably not much larger than that found in our 
surroundings, and is certainly very much smaller than that in NGC 
3115. 

The smallness of the dispersion in the radial velocities of nebulae 
in the Large Magellanic Cloud** indicates that the mass of this sys- 
tem can hardly exceed 2X 10°; since the luminosity is about 3X 10° 
times that of the sun, f would be 6 at most—again very much smaller 
than the values found in the outer regions of NGC 3115. 

No reliable data for other systems have yet been published. It 
should be noted, however, that measures by Horace Babcock*? of 
the rotational velocity in the outer regions of the great Andromeda 
nebula indicate that the average value of f for this system would be 
about 55. 

The dynamical models discussed do not give a direct explanation 
of the sharp edges of the isophotal surfaces of NGC 3115 near its 
equatorial plane. The most probable explanation would seem to be 
that these are due to the admixture of a relatively small fraction of 
stars with low velocities. 

In concluding this reconnaissance of the possible dynamical pe- 
culiarities of NGC 3115, it is well to emphasize its hypothetical char- 
acter. Of the hypotheses underlying the analysis the most uncer- 
tain is no doubt the neglect of absorption within the nebula and the 
supposition that the light we observe comes directly from stars. If 
important absorption and dispersion of light has occurred, this 
would probably have tended to lessen the light-gradients; the true 
disparity between the distribution of luminous stars and of mass 
would then be still greater than that indicated. The numerical re- 
sults for the velocity distributions would, however, lose their signifi- 
cance. The possibility that the equatorial plane of the nebula may 
be somewhat inclined to the line of sight also introduces some uncer- 
tainty; however, the sharpness of the edges of the isophotes near the 
equator shows that the inclination can hardly exceed 5°. The prin- 
cipal excuse for attempting an analysis in the present stage is that it 
may give some indication of what kind of additional observations 
would seem most promising for a better understanding of the struc- 

33 R. E. Wilson, Lick Obs. Pub., 13, 189, 1917. 

39 To be published in the Lick Obs. Bull. 
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ture of a system of this kind. There cannot be any doubt that an 
extension of the measures of rotation to greater distances from the 
nucleus would be of exceptional interest. Also, our knowledge of the 
system could certainly be greatly advanced if it were possible to 
obtain spectroscopic estimates of the true mean peculiar velocities. 


I am greatly indebted to Dr. Adams, who permitted the use of the 
60-inch telescope for the work on the light-distribution in nebulae. 
Likewise, my sincere thanks are due to Dr. Oosterhoff for his willing- 
ness to take the plates. To Mr. Humason I am indebted for showing 
me unpublished details of his measures of the rotation of NGC 3115. 
Professor Ornstein kindly permitted us to use the Moll registering 
photometer of the physical laboratory at Utrecht. Throughout, 
Messrs. Kleibrink and Pels have most efficiently co-operated in the 
measurements and computations. 

ASTRONOMICAL OBSERVATORY 


LEIDEN, HOLLAND 
October 23, 1939 
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IX. ON THE EXCITATION OF FRACTIONAL MULTIPLETS 
BY ELECTRON CAPTURE 


GEORGE H. SHORTLEY AND DONALD H. MENZEL 


ABSTRACT 

Observations show that certain permitted multiplets of O m1 and N 1m are not com- 
pletely represented in the spectra of gaseous nebulae. Some mechanism apparently acts 
to give selective excitation of certain levels of a term, without exciting other levels. 
Bowen attributes the phenomenon to resonance absorption, principally of the ultimate 
line of He 11, which happens to coincide with an ultimate line of O m1. The absorptions 
originate from the highest, instead of the lowest, levels of the ground terms. Since these 
levels have mean lifetimes of the order of several hours, Bowen’s theory requires the 
operation of some additional mechanism to maintain a sufficient population in the 
levels involved in the fluorescence process. Another method of excitation of fractional 
multiplets, direct electron capture by an ion in the ground level of a term (e.g., ?P1/2 of 
O Iv), is examined. The components predicted by quantum mechanics do not, however, 
agree with those observed in the nebulae. The conclusion is reached that Bowen is 
probably right, although a few discrepancies still remain unexplained. His method ap- 
parently requires inelastic electron impacts to be sufficiently frequent to maintain ap- 
proximately a thermodynamic equilibrium population of atoms in the levels of the 
ground term. This paper includes the relevant quantum mechanical analysis of the rela- 
tive probabilities of capture into the different levels of a configuration. 

Several years ago William H. Wright’ identified a number of lines 
in the far ultraviolet spectrum of the planetary nebulae with certain 
permitted lines of doubly ionized oxygen. A most striking feature 
of the identification was the fact that the complete multiplet 
2p3p *P — 2p3d 3P° does not appear, but only the lines arising from 
the level 3P°. 

In Table 1 are given the lines of two significant transition arrays 
of On. The lines observed in the gaseous nebulae are shown in 
boldface type. Nebular spectra to the violet of about \ 3200 are only 
incompletely observed. With the exception of the lines \ 3122 and 
d 5592, all of the observed lines would emanate from initial excita- 
tion of the level 2p3d*P%, with subsequent cascadings to lower 
levels. 

Bowen’ has suggested a possible mechanism for the selective ex- 
citation of this single atomic level. He has noted that a very close 
coincidence exists between the frequencies of the ultimate lines 
1s’*S — 2p?P° of Het and 2p? *P, — 2p3d 3P$ of O m1. Most of the 

* Pub. A.S.P., 46, 142, 149, 280, 1934. 2 Pub. A.S.P., 46, 146, 1934. 
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nebulae where fractional multiplets of O ut are observed are evi- | 
dently of high excitation. The great intensity of the line \ 4686 of 
He tt suggests that the ultimate line of He 1, referred to above, is 
also very intense. Bowen argues that resonance absorption by O 11 
of the intense He 11 radiation would lead to excitation of 2p3d 3P° 


TABLE 1 


TRANSITION ARRAYS 2p3d- 2P3P AND 2P3P- *2P3 s OF O III 
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& | 3p I jean 3774.01 | 3023. 42! 3035.92] 3043.01/ 3312.35 } 
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* Faint unobserved or unidentified line. Triplets according to ‘Mihul, Ann. de phys., 9, 261, 1928 


Singlets according to Edlén, Nov. Act. Reg. Soc. Sci. Upsala, IV, 9, No. 6, 1933. Edlén also identifies a num- 
ber of faint intercombination lines not included in this table 
t Although X 3122 does not originate in 2p3d sP2, Wright and Bowen have suggested that it may arise 4 
from resonance absorption of \ 303, by atoms moving with relative velocities of go km/sec f 
t Bowen and Wyse record the presence of \ 5592 in Nebulae NGC 7027 and 7662, with approximately 


one-half the intensity of \ 3760. We are indebted to these scientists for the communication of this informa 
tion in advance of publication. 


but not to excitation of the other levels of *P°, because the wave 
lengths do not coincide. As a further argument in support of his 
theory, Bowen points out that the line 2p3s *P? — 2p? 3P, of O m1, 
one of the end-products of the cycle of excitation and radiation, 
coincides in turn with the line 2p ?P°,, — 3d?D,,, of N ur. This co- 
incidence initiates another cycle, with further excitation of frac- 
tional multiplets of V 1. 

Bowen’s ingenious explanation, if true, places some definite re- 
strictions on the physical conditions obtaining in planetary nebulae. 
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It does not seem to have been emphasized previously that the levels 
3P, of Orn and ?P,,, of N ut, from which the original fluorescent 
absorption arises, are not the ground levels of the terms involved. 
Furthermore, although the mean lives of the levels 2p? P, and 
2p’P,,, are long,’ they are, nevertheless, much shorter than the 
expected time the atom must wait before being excited by the fore- 
going mechanism. Therefore, if only radiation processes are in- 
volved, we should expect the large majority of the atoms to be in 
the true ground levels, *P, and ?P,,,, respectively, of O mr and N 11. 
That a condition of this sort holds in interstellar space is shown by 
the observations of T. Dunham,‘ who points out that the only ob- 
servable interstellar absorption lines of 77 1 are those arising from 
the true ground level. Lines from other levels of the same term are 
not observed. 

In view of the fact that there is an alternative method by means of 
which fractional multiplets may be obtained—namely, in electron 
capture by ions not statistically distributed—we have thought it 
wise to undertake the relevant quantum mechanical analysis. Al- 
though the results are apparently negative, as far as explanation of 
the fractional multiplets of planetary nebulae is concerned, the calcu- 
lations appear to be of sufficient interest to warrant a summary here. 

The problem is extremely simple to visualize. For the sake of ex- 
plicitness, we shall discuss first the case of radiative electron capture 
by O rv to give directly O 11 in the configuration 1s? 2s? 2p3d. Under 
the conditions that obtain in the nebulae we might reasonably as- 
sume that practically all the O rv ions are in the very lowest level 
1s* 2s* 2p ?P,,,. If the O Iv ions were distributed between *P,,, and 
*P_,, according to the statistical weights, we should obtain (because 
of the sum rules) approximately a statistical distribution of the cap- 
tured electrons among the levels of O m1 2p3d. But on the assump- 
tion that all the ions are in the level ?P,,,, we expect and obtain, 
large peculiarities in the distribution in 2p3d. 

To make a (dipole) radiative transition into 2p3d, the electron 
before capture must be in either a p or an f orbit. The O Iv ion in 
the level 2p ?P,,,, plus the free electron in a hyperbolic orbit of energy 

3 The mean lives are: O m1 2p? 3P,, 2.9 hours; N 11 2p ?P;/2, 5.3 hours. 


4 Nature, 139, 246, 1937. 
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E, are to be considered together as forming a term of the continuous 
spectrum of O 1m, which we write as 


2P1/2EP1/2 , 2P1/2Ep3/2 , 2P1/2Efs/2 , 2pr/2Ef7/2 , (1) 


according to whether the free electron is in a p or an f orbit and ac- 
cording to the total (spin+orbital) angular momentum of this free 
electron. These are definite terms of O mr in jj coupling, capable 
of making radiative transitions to the levels of 2p3d, which are, how- 
ever, in LS coupling. The problem is then one of obtaining the 
strengths of the lines in the transitions from the configurations 2pEp 
and 2pEf in 7j coupling to the configuration 2p3d in LS coupling. 
This problem may be solved in all detail, but with considerable labor, 
by the standard methods of transformation theory; but a theorem 
we shall shortly prove will furnish very quickly the answers to the 
questions of interest here. 

We are interested only in the fofal strength from all the levels 
comprehended in (1)5 to a given level of 2p3d. /f the configuration 
2p3d were in jj coupling, the answer would be given very simply by 
the sum rule:° the total strength to any level (2p,,, 3d; J) in which 
the 2p electron has the same j value as in (1) would be proportional 
to 2] + 1; the strengths to all levels for which the 2p electron has 
the j value } would vanish. The theorem below shows that with 2p3d 
in any other coupling, in particular Russell-Saunders, we need only 
express the eigenfunction of each level in terms of 77 coupling 
levels, e.g.,7 





=) la V28 Mrady) V14¥(Pr/ads/2) 
ss 2V 21 W(py2dy/2) — 4¥(Psyads 2) , 
HpdP.) = (T= \[— Vs HPuady) + ab Psadye) + 30 Pvede)] 


¥(pd3F,) = (; 


etc. We then take the sum of the squares of the coefficients of the 
states p,,,d,,, and p,,,d,,, (which sum we call the p,,, parentage factor) 
$T.e., the levels 2p,/. Ep:/2, J = 0, 1; 2px/2 Eps/2, J = 1, 2; etc. 


6 Condon and Shortley, Theory of Atomic Spectra, p. 279. 


7 [bid. All the necessary transformation coefficients are tabulated on p. 294. 
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and multiply this sum by 2J + 1 to get the relative total strengths 
of the capture transitions to the different levels of 2p3d. For ex- 
ample, for *F, and 3P, the p,,, parentage factors are 


I 


(1/225 + 2)(16 + 21 + 14) = jf, 


(3'0)5 — b ’ 


respectively. The relative total strengths are then } + 5 = 35/9 


and} - 3 = 1/2. 

It is in this way that the p,,, parentage factors and the relative 
total strengths of Table 2 are computed for various transition arrays 
of interest.* These values are independent of the energy and the 
angular momentum of the captured electron, but to get transition 
probabilities from these strengths requires multiplication by the 
cube of the wave number of the radiation emitted in the capture, 
averaged in a complicated way over the energy of the captured elec- 
tron. The transition probability falls rapidly with increasing elec- 
tron energy. Consequently, for low electron temperatures, 7, such 
that kT < x, where x is the ionization energy, one should get satis- 
factory relative transition probabilities by using the cube of the 
wave-length limit of the continuous radiation. For very high tem- 
peratures the averaged factors for the various continua approach 
equality. 

The method used in these computations is a special case of a rather 
general theorem of matrix mechanics, which is demonstrated at the 
end of the paper. For the modifications necessary to handle capture 
into a configuration, such as p’, involving equivalent electrons, see 
the discussion under this theorem. We shall proceed first to the 
astrophysical discussion of Table 2. 

In Table 3 are given the relative strengths of the lines in the con- 
figuration arrays ps — pp and pp — pd. The relative intensities are 
to be obtained from these strengths by multiplication by the fourth 
power of the wave number and by a population factor that is unity 
if each upper level is populated in proportion to its statistical weight 

* If we had been interested in the case in which the ion is in ?P,,,, we should, by the 
same procedure, obtain the complementary relative strengths, which with the strengths 
from ?P,/, make up the whole 2J + 1. These are also listed in Table 2. 
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(2J] + 1). The population of a level depends on the rates at which ) 
atoms enter and leave the level, by all processes. For the assumption | 
of excitation of pp and pd, respectively, purely by electron capture, 

















TABLE 2 | 
RELATIVE TOTAL STRENGTHS OF TRANSITIONS INVOLVED IN RADIATIVE 
ELECTRON CAPTURE BY AN ION IN LEVELS p “P;/2 AND p?P;/2 
Ton in LEVEL 
TRANSITION TERM Sum Pe/s Paes 
AGE FACTOR 
*Pi/2 ?P3/a 
Capture of p 3F, ° 9 9 ° 
or f electron sF; 28/9 35/9 7 4/9 
to give pd 3F, 35/9 10/9 5 7/9 
'F; 7/3 14/3 7 1/3 
3D; 14/9 49/9 7 2/9 
3—, 35/18 55/18 5 7/18 
3D, 3/2 3/2 3 1/2 
'D, 5/3 10/3 5 1/3 
3P, 5/2 5/2 5 1/2 
IP, 1/2 5/2 3 1/6 
3P, ° I I ° 
tp, I 2 3 1/3 
Capture of s iD, ° 7 7 ° 
or d electron 3D), 5/2 5/2 5 1/2 
to give pp 3D, 5/2 1/2 3 | 5/6 
'D, 5/3 10/3 5 1/3 
iP, 5/6 25/6 5 1/6 
3P; 3/2 3/2 3 1/2 
iP, 2/3 oy I 2/3 
i¢ I 2 3 1/3 
3S; I 2 3 1/3 
™So 1/3 2/3 I 1/3 
Capture of s 3p, 5/3 25/3 10 1/6 
or d electron 3P, 3 3 6 1/2 
to give p? 3P, 4/3 2/3 2 2/3 
‘1D, 10/3 20/3 10 1/3 
So 2/3 4/3 2 1/3 
Capture of p sP, fe) 5 5 fe) 
electron to 3P, 2 I 3 2/3 
give ps 3P, I ° I I 
'P, I 2 3 1/3 
| 

















Table 4 takes into account solely the rate of entrance of the atoms. 
It is obtained from Table 3 by multiplication of the files referring to 
the upper levels by the p,,/, parentage factors of Table 2.° 

9 Note also that the variation in the (wave number); factor necessary to convert 
Table 2 into transition probabilities is neglected. 
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Table 4 contains only part of the story, however. Since 2p3d feeds 
very strongly into 2p? as well as into 2p3p, only a fraction, and prob- 
ably a small fraction at that, of the atoms entering 2p3d will be 
available for transition to 2p3p. Most of the atoms in the terms 
3D, 3P, *F, *D, or 'P of 2p3d will have a much higher probability, 


TABLE 3 


LINE STRENGTHS IN NORMAL ARRAYS 



























































IS? 2s? 2p3p 

iD 3p 3$ *D Pp 1S 

J 3 2 I 2 I I 2 I ° 

4 54 ° ° ° fe) ° ° ° ° ° 

3F 31 14 112/3 ° ° ° fe) fe) ° ° ° 

2 2/15 14/3 120/5 re) fo) ° ° ° ° ° 

3 | 28/3 7/6 ° 63/2 ° ° ° ° ° 

3—p 2 7/6 | 125/24 9/8 45/8] 135/8] o ° ° ° ° 

A I fo) 9/8 27/8 3/8] 45/8) 15/2| o ° ° ° 
a f | ses cai a ee ee 
h 2| 7/10} 1/8 1/120] 75/8] 25/8] o | 50/3| o ° ° 

N > > 
, i I fe) 3/8 1/8 25/8} 15/8] 5/2| 10 ° ° fe) 
- ° ° ° 1/6 ° 5/2] oO 10/3} oO ° ° 
3 | Gis tole ae a 
if | 3 | fe) re) ° fe) | ° ° ° 42 ° ° 
: Me Se See: 
‘DPD | 2 | ° ° | o | o | ° ° ° 15/2| 45/2] o 
| — — ee | a _ <2 
iP | I | ° | o | oO | ° ° | ° 1/2} 15/2] 10 
2-4 aa We ee, aly a 
a 2 7 5/4 1/12 | 15/41 5/4] 0 5/3) oO ° ° 
a | 3P I ° 15/4 ela Ps 4| 3/4) 1 I ° ° ° 
, ° fo) ° 5/3 | ° | I ° 1/3| oO ° ° 
: 3 ess SR, ON! bars 
ms iP | 1 | ° | ° | ° ° | ° ° | ° | 5 3 I 
| 








because of the large v-cube factor, of relapsing directly to 3P, "D, or 
'S of 2p’ than of going to the terms of 2p3p. Atoms in any one of 
the levels of *F of 2p3d, however, can relapse only to *D of 2pa3p. 
The intensity of this multiplet remains as given in Table 4, whereas 
all the others are diminished by a very appreciable factor. The 
strengths in the second part of Table 4, from the array 2p3p — 2p3s, 
are unaffected by any such complications. Comparison of Tables 3 
and 4, even without the depletion just discussed, shows that the 
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relative strengths of the various multiplet lines are abnormal. Cer- 
tain lines are missing or are very weak. As far as the multiplets 
3P — 3[—, 3P — 3P, and 3P — 3S are concerned, the intensity distor- 


Te 


TABLE 4 
MODIFICATION OF TABLE 3 ON THE ASSUMPTION OF EXCITATION BY DIRECT 
ELECTRON CAPTURE IN THE UPPER LEVELS OF 2p3d—2p3p AND 2p3p—2p3s 

































































‘D iP S| Dp | P| ss 
J | 
3 2 | I 2 I ° I 2 I ° 
Is? 28? 2p3p 
> | 
af 4 ° ° ° ° ° fe) ° ° ° ° 
2 3] 3 | 56/27 | 448/27 ° ° ° ° Ot O11 O11 °o 
a 2 | 14/135} 98/27 98/5 ° ° ° ° ° ° ° 
2) 
> cs es 2 
a 3 | 56/27 | 7/27 | ° 7 | ° ° ° ° ° ° 
z 3D 2 | 49/108) 875/432] 7/16 | 35/16) 105/16} o ° ° ° ° 
= I o | 9/16 27/16 3/16] 45/16]15/4 | o| of of Oo ) 
i, ee ee en: ED ESAS eta eee 
= | Ree 
= 2 7/20 1/16 | 1/240) 75/16) 25/16] © [25/3] © ° ° 
™ | 3P I ° 1/16 | 1/48 | 25/48} 5/16] 5/12| 5/3} oO | o| o 
= ° eS ° . T 2 | o | 0 ° | o| o| o 
%, et ae Se 4 = oh, eer eee. eos | 
nN | | | 
% } 3 fr) | fe) | fo) | o | 0 | | ° | 14| o| o 
1D 2 ° | ° | ° | | ° | ° | ° | 5/2|15/2| o 
ay a B aA . | c ie: SY ere 
I I ° ° | ° | ° | ° ° Oo} 1 . 5 : me 3 
} 
1s? 2s? 2p3p (excited by capture) ) 
2 
4 | 
% 2 ° | 5/8 | 5/72 5/8 5/8 ° 5/9} oO ° ° 
% 3P I o | 15/8 | 25/24 5/24 3/8 | 2/3) 1 3| ° ° ° 
ns ° o | o | 25/18 ° 1/2 ° I 0| ° ° ° 
ip | I Page 4 ° o | o ° | ° ° | 5/3 ae ie 
| | | | 5 | 3 

















tion agrees fairly well with observation. The partial excitation of *P, 
is accounted for naturally, whereas Bowen was forced to consider the 
extremely high velocity differential of go km/sec to achieve excita- 
tion of 2p3d %P, by resonance absorption of the ultimate line of 
He u. But failure of at least the lines of the *F — *D and *D — 3P 
multiplets to appear is definite evidence against the validity of the 








en 
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mode here considered of exciting fractional multiplets. It is true 
that we have taken into account only the excitation of 2p3d directly 
by electron capture, while electrons will also be captured into higher 
configurations and feed down into 2p3d. But the changes occasioned 
by taking this into account will certainly decrease the intensity 
anomalies and will not aid matters. 

Looking now at 2p3p — 2p3sin Table 4, we see that direct capture 
into 2p3p could not possibly account for the observed intensity 
anomalies. A rough computation shows that the assumption of ex- 


TABLE 5 
RELATIVE STRENGTHS OF OBSERVABLE LINES ON BOWEN’S HYPOTHESIS 












































Is? 2s? 2p3p 
J ‘D sp 3S 
| 

3 2 I 2 I ° I 
1S? 2s? 2p3d. | 3p | 2 7/10 | 1/8 | 1/120 75/8 25/8 ° 50/3 
2 7/10 | 1/32 ly ‘4320 | 225/32 | 125/96 | Oo 250/27 
1s? 2s? 2p3s..| 3P I fe) 3/32 | 1/288 75/32 75/96 ° 50 9 
° ° ° 1/210 ° 25/24 ° 50/27 





citation of 2p3p principally by transitions from 2p3d does not im- 
prove the agreement. 

Also, the observed anomalous intensities in the multiplet 3p ?P° — 
3d?D of N m1 cannot be accounted for by the present scheme, be- 
cause the ground term of JN Iv is 'S,, a single level. 

We are apparently forced to return to the Bowen mechanism for 
the excitation of these lines, despite the fact that it does not account 
for the strength of the singlet line \ 5592. In Table 5 are given the 
relative theoretical strengths of the lines of the significant multiplets 
of O1n, calculated on the basis of initial selective excitation of 3P, 
alone. The strengths in the two transition arrays of this table have 
the correct relative values and need only be multiplied by v4 to give 
the relative intensities of all the lines. The tabulated strengths show 
good qualitative agreement with the lines observed in the nebulae. 

As was noted earlier in this paper, acceptance of the Bowen mech- 
anism places a definite restriction on physical conditions within a 
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nebula. In interstellar space the high concentration of atoms occurs 
in the lowest level, whereas in the nebulae an appreciable population 
must be maintained in the levels *P, and *P°,, of O m1 and N m1, re- 
spectively. The mean lifetimes of these levels are of the order of 
several hours,’ and a concentration will occur only when the rate of 
excitation to these levels is also about once in three hours. Radiation 
processes, with the high dilution existing in nebulae, cannot be called 
on for so frequent excitations. A study of the problem indicates that 
collisional excitations from the ground level must be sufficiently 
frequent to maintain the required population. This investigation, 
by Hebb and Menzel, will form the subject of the next paper of the 
present series. 

We proceed now to the proof of a theorem concerning the trans- 
formation properties of the line-strength matrix which will justify 
the method used in the above computations. 

THEOREM.—Consider a configuration a and all the levels BJ of a 
particular angular momentum J belonging to a configuration B. In a 
particular coupling let the levels of a and of BJ break into two groups 
a’, a” and BJ, B'S such that B’’'J does not combine with a’ (B’J may 
combine with both a’ and a’). In this case the matrix’? of S'/* connect- 
ing the levels of configuration a with those of BJ will have the form 











, ” 
B’J BJ 
| 
ae eee Sup ° Oo 
, Sar S22 - « Say 1) 
a 
Snt Sn2---. Sn» Oo O ce) 
(2) 
X X X Xx X X 
o a ee : x x 
a 
X Be te 58. a eae | X , eee 3 











t0 See Condon and Shortley, of. cit., p. 277, for definition of this matrix, whose com- 
ponents are the square roots of the line strengths chosen with phases corresponding to 
those of the matrix components of electric moment. 
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Let the strength of each of the files connecting one of the levels of BJ 
with all the levels of a’ have the same value Cy, and let all pairs of these 
files be orthogonal. That is, let 


al 
> se Sin = Cy ben . (3) 


1 I 


Now consider a level BJ of configuration B in any other coupling. When 
this level is expressed in the B’, B’’ scheme, let the sum of the squares of 
the transformation amplitudes connecting this level with the levels of 
B’J be tr. Then the total strength of all transitions connecting the level 
BJ and the levels of the part-configuration a’ will be rCy. 

To prove this statement let the level BJ be connected with those 
of B’J and B’’J by transformation amplitudes contained in the one- 
column matrix 











BJ 
t, 
B’J t 
bi 
—*— (4) 
x 
X 
p’"J 
x 
where by hypothesis 
= @=7r<. (5) 


Now, since a’ does not combine with 8’ J, we may get the S”? ma- 
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trix connecting a’ and BJ by multiplying the upper left square of (2) 
with the upper part of (4) to get 














The sum of the squares of these elements is 


n v v 
7 =m ~ 

SN S Sight Sint) = Cy ¥ #2 = 7Cy, 
yy iy a_i * 

t=1 =I n=1 g§=1 


where the right-hand — are obtained by carrying out first the 
sum over 7, then that over &. Note that in cases where this theorem 
is valid, one may obtain the desired sums without knowing the rela- 
tion between the phases in (2) and the phases in (4). This is some- 
times a great convenience. 

We may now apply this theorem explicitly to transitions from 
aE] to an’l’ where the electron n’l’ is not equivalent to any in the 
parent configuration a. Let a, represent a given /evel of the parent 
configuration in the ion. Then, if we characterize the levels of an’l’ 
by the same parents, there will be no transitions from any level of 
a:n'l’ to any levels of aEl except those of a,El. The sum of the 
strengths from a given level a,n’l’J to all levels of a,El is 


(2J] + 1)| (El: Pin’l’)|? (I, 1’) , (6)! 


11 [bid., p. 281. Here =(1,1+ 1) = (1+ 1)(21+ 3); Z(1,1 — 1) = I(2l — 1); (El: Pin’) 
is an integral of the product of 7 and the radial functions of the two states as defined in 
ibid., p. 132. Since one of the levels is in the continuum, we must consider (6) as re- 
ferring to unit energy range and the radial function for the continuum normalized in the 
proper fashion. See Kemble, Quantum Mechanics, for a good discussion of this normal- 


ization. 
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proportional to 2J + 1. The orthogonality conditions of the fore- 
going theorem are also satisfied; hence, in the case of a level of 
an’l’J in any other coupling, having a total component 7 in the sense 
of (5) in the direction of a,n’l’J, the total strength to all levels of 
a,E] will be 7Cy. Note that the relative values of the strength sums 
from the different levels of a given discrete configuration an’l’ to con- 
tinuous states a,El is independent of whether | = 1’ + 1 or I’ — 1, 
and of the energy E. 

O 11 is just a special case of the previous paragraph, with a = 
IS? 2s? 2p; a; = 1s? 2s? 2p *P, ,,; hence our earlier calculations are justi- 
fied except for the case p(?P,,,)Es or p(?P,/.) Ed > p’, which violates 
the inequivalence specification. But in this case an examination of the 
S'/? matrices in 77 coupling shows them to be of the form (2), if for the 
p’ configuration we take B’o = px2p1/20; B’I = P3/2P1/21; B’2 = Ps/2Pr/22 
(for the case ps, p? see the matrix on page 265 of Condon and 
Shortley). The Cy’s in these cases are no longer proportional to 
2J + 1 but have the same relative values for capture of s and d elec- 
trons, namely, Co:Ci:C, = 2:3:5. 
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AN INVESTIGATION OF THE ROWLAND 
INTENSITY SCALE 
DONALD H. MENZEL, LEO GOLDBERG, AND EDWARD M. COOK 


ABSTRACT 

A new calibration of Rowland’s scale of solar intensities has been derived with the 
aid of the theoretical strengths of multiplets in transition arrays. The calibration takes 
the form of a double-entry table giving average values of log Xo, where X, is the optical 
depth at the center of an absorption line, for each Rowland intensity from —3 to +8 
at intervals of 400 A from \ 2800 to \ 6800. A total of 37 transition arrays, representing 
13 metallic elements and 111g lines in all, was employed in the analysis. The present 
calibration is considerably less steep than the earlier one of Russell, Adams, and Moore. 

Probably the most promising method of deducing the physical 
conditions in the solar reversing layer is through studies involving 
the curve of growth of absorption lines in the Fraunhofer spectrum. 
Such studies have been greatly accelerated by the wealth of equiva- 
lent widths published by C. W. Allen,’ and further important results 
are to be expected on the completion of the Utrecht atlas of the solar 
spectrum. In the future we may expect these photometric measures 
to replace completely the admittedly rough estimates of intensity 
contained in the Revised Rowland Tables. Nevertheless, as ragged 
and approximate as the Rowland intensity estimates are, they repre- 
sent a vast amount of data. The estimates are particularly impor- 
tant for the very faint solar lines and for lines in the region 3000 
4000 A, for which measures by Allen are not available. 

Before the Rowland estimates can be employed in analyses of the 
physical conditions in the solar reversing layer, it is necessary that 
the arbitrary intensity scale be calibrated, preferably in terms of a 
quantity related to the numbers of atoms per square centimeter 
above the photosphere that are active in producing the various 
Fraunhofer lines. Such a calibration was first carried out by Russell, 
Adams, and Moore’ in their pioneering investigation of a decade ago. 
The important astrophysical results that have accrued from this 
initial investigation are so well known that there is no necessity of 
reviewing them here. 

The Russell-Adams-Moore calibration assumed that the relative 

* Mem. Comm. Solar Obs., 1, No. 5, 1934. 7Ap. J., 68, 1, 1929. 
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strengths of the lines in a multiplet correspond to the relative num- 
bers of atoms that are active in producing these lines. For each mul- 
tiplet the lines were combined into two or three groups, for which 
mean values of log N and R yielded one or two values of AR and 
A log N corresponding to a given value of R. The limitations of such 
a calibration are obvious. Confining the theoretical strengths to lines 
of the same multiplet meant that only a relatively small range of in- 
tensity could be covered in any one group of related lines, a pro- 
cedure that increased the possibilities of systematic and accidental 
errors in the calibration. Also, since lines of the same multiplet do 
not differ widely in wave length, there is no possibility of determin- 
ing whether a line in the red requires more, or fewer, atoms for its 
production than one of the same Rowland intensity in the violet. 
We may overcome these two difficulties simultaneously by em- 
ploying the theoretical strengths of different multiplets in transition 
arrays. Whereas a single multiplet rarely consists of more than 12 
or 13 lines with approximately a 100 A wave-length range, the num- 
ber of lines that can be related in a transition array may number well 
over a hundred and possess a wave-length dispersion of several thou- 
sand angstroms. The relative strengths of different multiplets in a 
transition array are calculable according to methods outlined by 
Kroénig,? Shortley,4 Condon and Ufford,’ and Menzel and Goldberg,° 
under the assumptions that configuration interaction is negligible 
and that the atom obeys the rules of LS coupling. These idealized 
conditions are not too well fulfilled by complex atoms, and the pres- 
ent investigation was begun with the hope of testing the general 
validity of the theoretical multiplet strengths. But with the progress 
of this and other related solar work,’ the evidence accumulated that 
the use of large numbers of multiplets tends to iron out random 
deviations and yields consistent results in the mean, although the 
theoretical strength of a single multiplet may be considerably in error. 
Consequently, the decision was made to employ the theoretical 
multiplet strengths in an independent calibration of the Rowland 


scale. 
i Zs. f. Phys., 33, 261, 1925. 5 Phys. Rev., 44, 740, 1933. 
1 Phys. Rev., 47, 295, 1935. ° Ad. J., G4, 1, 1936. 


’ Menzel, Baker, and Goldberg, Ap. J., 87, 81, 1938. 
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Menzel* has shown that the quantity W/X, where W is the equiva- 
lent width of an absorption line and } its wave length, is a simple 
function of X,, the optical depth at the center of the line. The ex- 
pression for X, is written 

Xo = 3.510 X 10°” ae eX yf kT e 4s =. (1) 

: b(T) ro SS" 
where JN, is the total number of atoms of a given element in a par- 
ticular stage of ionization in a square-centimeter column above the 
photosphere, b(7) is a partition function, k is Boltzmann’s constant, 
x, is the lower excitation potential, 7 is the temperature, u is the 
atomic weight, and s/Zs is the fraction that a line in a given multi- 
plet contributes to the total strength, S, of the multiplet. The 
strength is measured in units of 0”, the square of the first-order elec- 
tric moment between the two electronic configurations. ¢ is a di- 
mensionless constant, defined by the relation 
0 


RIO: er (2) 


ae 





where d, is the radius of the first Bohr orbit, and ¢ is the electronic 
charge. 

In the interest of consistency with equivalent-width investiga- 
tions, we have chosen to calibrate the Rowland scale in terms of Xo. 
The quantity @ is constant for a given transition array. Hence, if 
we set 


log X%’ = log S — log = = = xy + 3 log xu, (3) 





X%’ proves to be proportional to X, for all the lines of a transition 
array. The quantity 7 is the excitation temperature governing the 
distribution of atoms among the various excited states. In a pre- 
liminary calibration, the value T = 5740° was adopted, but later 
investigation? showed that a temperature of 4500° corresponded 


8 Ap. J., 84, 462, 1936. 
9 Menzel, Baker, and Goldberg, op. cit.; R. B. King, Ap. J., 87, 40, 1938. 
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more closely to the observations, and this value was employed in 
the final calibration. It will be seen later that the slope of the cali- 
bration-curve is not very sensitive to the temperature assumed. 

The procedure followed in effecting the calibration may be briefly 
outlined. For all the lines of a transition array we calculate log X,’ 
from equation (3). The quantities S and s/s are taken from Gold- 
berg’s’® and from Russell’s" tables, respectively. A total of 37 transi- 
tion arrays, representing 13 metallic elements and 1119 lines in all, 
was employed in the analysis. The material is contained in Table 1, 
where successive columns show the element, the transition array, 
and the numbers of multiplets and lines in each array. Since the 
lines of a transition array usually extend over a considerable re- 
gion of the spectrum, we have attempted to establish the variation 
of log X;’ with \ for each Rowland intensity. The spectrum was 
divided into successive wave-length intervals and the average taken 
of all the values of log X%’ for lines of the same R that fall within a 
particular interval. It is, of course, desirable that the wave-length 
interval be as small as possible, so that rapid variations in Rowland’s 
scale may be detected. In practice, however, we are limited by the 
relative thinness of the distribution with wave length of the set of 
related lines. Accordingly, we found it most practical to choose in- 
tervals of 400 A over the region from \ 2800 to \ 6800. The cali- 
bration for each array then takes the form of a double-entry table 
which gives a value of log X4’ corresponding to each value of R and X. 
Table 2 provides an illustration of the procedure, showing the two- 
dimensional table for the transition array Ti 1 3d°4s — 3d%4p. 

If a sufficient number of lines in each array were available, each 
such table would represent a two-dimensional calibration of Row- 
land’s scale. It would then be a simple matter to obtain the final 
calibration by combining the various transition arrays graphically. 
Unfortunately, however, although the maximum possible number 
of entries per table is 144, the actual numbers vary from as few as 
2 or 3 to only 38 for the transition array shown in Table 2. We must 
therefore seek to combine the various tables by some numerical 


10 Ap. J., 82, 1, 1935; ibid., 84, 11, 1936. 


1 Ap. J., 83, 129, 1936. 
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TABLE 1 


SUMMARY OF SPECTROSCOPIC DATA 








Element 








Cat. 


oct 


Sc II 


tS: re 








mae Number of 
Transition Array Multiplets 
4s4p — 484d I 
484p — 4p” 3 
484p — 486s 2 
484p — 486d | 2 
484P — 4878 2 
4s3d—4s5p | 2 
3d4s— 3d4p 4 
483d —4s6p | 2 
3d4s?— 3d4s4p | 2 
3d24s — 3d24p 6 
3d4s—3d4p 6 
3d?— 3d4p 8 
334s — 3d34p 34 
3d74s?— 3d?4s4p 18 
3034s — 3d74s4p 9 
3d24s4p — 3074858 8 
3d74s4p — 3d?4s4d II 
3d3— 3d24p 24 
3d?4s — 3d74p 17 
3d44s— 3d44p II 
3d54s— 3d54p 2 
3d*4sgp — 3d 4485s 5 
3d44s— 3d ‘4p 3 
3d°4s— 3d°4p 6 
3d54s4p — 3054858 3 
3d54s4p — 3d54s4d 2 
3d74s— 3d74p II 
3d°4s4p — 3d°4s5s 8 
3d°4s— 3d°4p 5 
3d*4s4p — 3d*4s4d 13 
3d°4p — 3d°4d 16 
3d°4p — 30%6s 4 
3d%4p — 3d%5d 4 
4d?5s?— 4d?5s5p 6 
4d35s—4d35p 5 
4d?5s—4d?5p 14 
4d3—4d?5p 14 











Number of 
Lines 


6 


oo) 


NS 
mann fw 


54 
II 
13 


64 
82 
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method.” The following method of effecting the combination ap- 
peared to be most satisfactory. 

Let us denote an entry in any double-entry table by F’(R:A;),, 
where R; and ); are particular values of R and X, and n refers to a 
given transition array or double-entry table. The corresponding 


TABLE 2 
MEAN VALUES OF LOG X@’ FOR TRANSITION ARRAY 77 I 3d34s—3d34p 














Wave LENGTH 

R 
4000-4400 | 4400-4800 | 4800-5200] 5200-5600] 5600-6000 | 6000-6400 | 6400-6800 
“4... .|-o. 24:€2)" Fo: 32 Ca) [-Fora6 GE: «soc e de Shee 5 osc opto ce eaeles wee es 
+3 —o.58 (1) |—0.45 (6) |—0.04 (4)].........].....006. ee) | 
+2....}—0.61 (7) |—0.67 (11)}—0. 58 (6)| —0. 62 (3)]......... — 6.40 Ca)E eo. 6 ces 
+-3r,...[—21.92 63) |= 8: tao) [8 OGG. 5 i crc checcvedevcbeaccovssd —o.35 (1) 
o..../—1.61 (7) |—1.71 (7) |—1.38 (4)... ..|—1.24 (1)|—1.15 (1)}—0. 50 (1) 
—1..../—1.75 (9) |—2.07 (8) |—1.69 (5)}—1. 88 (3)|......... cat Ce ft | eee 
ase TY Cee ee —2.27 (2) |—2.19 (6)/—1.99 (3)]......... et Ot) | rer 
a es) eee Ay mie SO Ca) b—2. EBC ns cis a Se ev eae apewacceuns —1.52 (1) 
































* Numbers in parentheses denote number of lines comprising each average. 


value of log Y, may be written F(R;A,;), since log XY, depends only 
on R and \. Equations (1) and (3) show that F’(R;\,;), differs from 


F(R;d;) by a constant, which we denote by the symbol A,: 





As = F’(Rid;) n a F (Rid) ) 
x =e Na (4) 
= 13.35 log ES 6]. : 
or 
F(RiAj) = F’(RiAj)n — An- (5) 


If we write equation (5) for any other transition array m, inserting 
the same values of R and X, we have 


F(R3X;) pa F’(Rid;)m was Be . 


(6) 


2 We are indebted to Dr. T. FE. Sterne for suggesting the proper method of approach 


to this problem. 
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Subtracting equation (5) from equation (6), we find 
F'(RAj)a — F(RA)a = On — Mn (7) 


Thus, any entry in one table theoretically differs from the corre- 
sponding entry in another table by a constant, and, if a single value 
of A is known, the remaining values may be calculated from equa- 
tion (7). Because of errors in the theoretical strengths, however, and 
also in the Rowland estimates, the differences between correspond- 
ing entries in any pair of tables are not identical. Some method of 
averages must be employed. To calculate the A,’s we begin by as- 
signing an arbitrary value of A, to array No. 1, 7i 11 3d° — 3d’4p, 
and then proceed to determine relative values of A,. Let 


Amn = Am si An ’ (8) 
whereupon equation (7) becomes, for m = 1, 
F'(R;): — F’(RA;)n = Am- (0) 


We can calculate a value of A,, for every pair of values F; and F/ that 
occurs in the same R, X interval. Since they are not identical, the 
simplest procedure, which immediately suggests itself, is first to 
compute the average value of A,, for each transition array and then 
combine the 37 tables by a straight averaging process. Unfortu- 
nately this method fails because a great many tables have few or no 
entries in common with the first table. Hence, we must first endeav- 
or to enlarge this table more nearly to completion. This result may 
be accomplished by combining the 6 most complete tables (Tables 
1-6) by a least-squares solution. 

We are required to evaluate the quantities A,., A,;,...., Aw by 
the method of least squares. If A,,, is the mean value of A,,, as com- 
puted from pairs of tabular entries with the same values of R and X, 
the A’s are to be computed according to the following equations 
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of condition, weighted according to the square root of the number 
of values of A,,, that were employed in computing the means: 





An = Mn (14) | 
ers (Via) | 
Au = Ang (Vs) | 
Ais = Ais (v8) | 
Aw = Ais (W'16) | 
An = —AutS,= hb, (Via) | 
An = —AntAy=Sy (V3) 
As = —AntAy = Az (V11) (10) 
Aw = —AntdAwe= As  (Wro) | 
Ay = — Ay + Ay = Ay (V8) | 
Ass = —Arzz+ dis = A;; (Vrs) | 
Ay = —A;+tAw=As (Vo) | 
As=—Autds=Ag (V4) | 
Ay = — Aut Aw = Ay (V3) | 
Ag = —Astds=ds (Vo) | 


parentheses. The equations were 
, Aw. The results are shown 


The weighting factors are given in 


solved in the usual way for A,., Ay;,... . 
TABLE 3 


RESULTS OF LEAST-SQUARES SOLUTION FOR Ayn 




















Table | Element Numer 
| of Lines 
I | Ti 1 3d3—3d24p 72 
2 Tit 3d?4s?— 3d74s4p 57 
3 Ti i 3d34s— 3d34 128 
|} 4813 p 
4 Fet 3d°4s4p—3d°4s5s 82 
or | V 1 3d44s—3d44p 60 
6. | Mnt 3d°4s—3d°4p 54 





in Table 3 





Arn Ain 
0.00 0.00 
—0o.96 —<0.07 
—0.04 —©. 160 
+2.81 +2.93 
2500 —1.02 
+1.73 +1.75 





Ie H+ H+ + 


0.0096 





. The second column contains the element and transition 


array; the fourth column, the quantities A,,; and the fifth and sixth 
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columns the computed values of A,, with their mean errors. The 
agreement of the A’s with the mean values and the smallness of the 
mean error speak well for the consistency of the results. 

The values of F’(R;\;), were next increased by the addition of the 
value of A,, to reduce the other five tables to the scale of the first. 
Next, all entries of Tables 1-6 with the same values of R and X were 
averaged. We then incorporated the remaining transition arrays into 
the combined table by computing the mean value of A,, for each 
table, adding this quantity to the entries F’(R;A,;), and averaging as 
before. The final calibration is given in Table 4 and represents the 
material from 37 transition arrays. The table gives a value of 


log Xq’ = log X, + A, , (11) 


where the subscript 1 refers to the array 77 11 3d’ — 3d’4p, for every 
pair of values of R and X, from R = —3 to R = +8 and from 
A 2800 to dA 6800. 

The calibration contained in Table 4 involves the assumption that 
the quantities A, are constant with wave length and hence implies 
that \,, the total number of atoms above the photosphere, is inde- 
pendent of wave length. If the general opacity of the solar atmos- 
phere varies with wave length, the position of the effective photo- 
sphere will rise or fall, and the numbers of atoms will also vary. To 
test for this variation and at the same time to impose a general check 
on the calculations, we have employed Allen’s measures of equivalent 
widths in an independent calibration for the region AX 4000-6600. 
Values of log X, were obtained for each unblended line in Allen’s 
“Intensity Table” from the tabulated equivalent widths and the 
curve of growth. We then computed mean values of log XY, for each 
Rowland intensity and for wave-length intervals of 400 A, to make 
the results correspond with the multiplet calibration. The results 
are given in Table 5. We examined the constancy of A, as a function 
of wave length by computing mean values of this quantity from 
Tables 4 and 5. The resulting values are plotted against the wave 
length in Figure 1. The straight line, which was fitted to the plotted 
points by least squares, obeys the relation 


A, = 2.72 — 0.230A, (12) 




















THE ROWLAND INTENSITY SCALE 329 


where } is in units of tooo A. This law corresponds very nearly to 
a \3 opacity variation. It must be remembered, of course, that the 
relation (12) may be due in part to systematic wave-length errors 
between the two calibrations. 


TABLE 4 


MEAN VALUES OF LOG X@ AS A FUNCTION OF R AND A 












































Wave LENGTH 
R 
2800-3200 3200-3600 3600-4000 4000-4400 4400-4800 
So eee ROE Nat +1.90 (1)*| +1.22 (4) | +1.19 (2) | +0.77 (3) |........... 
7 eataewaes +2.06 (1) | +1.53 (4) | +1.24 (2) | +0.45 (2) | +0.50 (1) 
PER 5 Scr pean +1.51 (1) | +0.92 (5) | +0.83 (2) | +0.88 (3) | +0.53 (2) 
“ptr Pee oe +o0.44 (1) | +0.48 (8) | +0.21 (4) | +0.03 (7) | +0.49 (3) 
i ene +0.87 (2) | +0.46 (8) | +0.19 (7) | +0.15 (8) | —o.04 (4) 
+3 -.e.s.s| +0.45 (3) | +0.06 (9) | —0.56 (6) | —0.40 (8) | —0.51 (0) 
+2 ....| —0.09 (4) | +0.06 (1o)| —o.69 (10)} —0.84 (9) | —o.91 (10) 
+1 .....-| —0.45 (3) | —0.42 (10)} —1.05 (11)} —1.11 (11)} —o.82 (8) 
O° -.eesee] —0.69 (5) | —0.90(5) | —1.51 (5) | —1.67 (8) | —1.23 (11) 
a ee ..seee{ 4.33 (4) | 1-42 (5) | —1.34 (4) | —2.27 (6) | —1.69 (9) 
—2..... .. ee] —1.28 (2) | —1.33 (4) | —1.42 (4) | —2.13 (5) —g 82 (6) 
at SEs | eerie: —1.64 (5) | —1.42 (3) | —2.15 (1) | —2.25 (3) 
Wave LENGTH 
R eioniceibinie —— 
4800-5 200 5 200-5000 5600- 6000 6000-6400 6400-6800 
+8 +0.89 (1) | , et tS eee +0. 37 (1) 
+7 +o0.74 (1) | +0.46 (1) Se ee rs) | Oe ee 
+6 +o.64 (1) | +0.18(2) | —o.15 (1) | —0.36 (2) | +0.04 (1) 
+5 +o.51 (2) | —o.21 (2) | —0.72 (1) | —0.69 (1) | —0. 26 (1) 
+4 +o.11 (3) | —0.59 (3) | —0.55 (1) | —0.82 (1) | —0.68 (1) 
+3 —0.52 (7) | —0.32 (4) | —0.69 (1) | —o.62 (2) |.....-..... 
+2 —o.73 (6) | —o.78 (4) | —3.35 (7) | ~0.86(%) I... 2.4.2. 
+1 —1.11 (4) | —1.56 (2) | —1.38 (4) | —1.55 (3) | —1-03 (3) 
o | —1.43 (7) | —1.36 (6) | —1.85 (8) | —1.59 (3) | —1.40 (2) 
—] ab 52 (6) | —1.72(7) | —1.88(6) | —2.09(7) | —2.31 (1) 
—2 | —2.00 (4) | —1.75 (8) | —2.01 (5) | —1.94 (5 —2.70 (1) 
oat | —1.99 (6) | —2.19 (6) | —2.97 (3) | —1.15 (1) | —2.49 (2) 
| | | 














Possible systematic differences between the two calibrations as a 
function of Rowland intensity should be revealed in the averages of 
log X, and log X%’ for each Rowland intensity from —1 to +6 over 
the wave-length range from 4000 A to 6400 A. The two sets of values 
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TABLE 5 


VARIATION OF LOG X,, OBTAINED FROM EQUIVALENT WIDTHS AND 
THE CURVE OF GROWTH, AS A FUNCTION OF R AND A 















































Wave LENGTH 
R = 
4000-4400 | 4400-4800] 4800-5200] 5200-5600] 5600-6000 | 6000-6400 | 6400-6800 
+6...... +2.80 2.09 +2.26 +2.10 +1.38 +0.87 +1.37 
5 +2.23 +1.84 +1.88 +2.04 +0.76 +0.70 +1.10 
+4 +1.82 +1.49 +1.79 +1.61 +0.73 +0.52 +o.82 
+3 +1.58 | +1.22 +1.21 +1.16 +o0.42 +0.36 +0.53 
a, +1.02 | +0.86 +0.85 +0.72 +o.18 +0.13 +0.27 
TPiecsix| Or02 +0.44 +0.39 +o. 28 —0.07 —O.14 0.00 
O......| +o.08 | —o.10 +0.03 —0.15 —0.39 —0.45 —0.44 
ee ai | “0.42 | =—O.50| =O. 3% =6. 50 => OF —©.00 =O. 52 
| 
20 | | | | | | 
° 
* 
IS - " ° = 
e o 
10 i | l ! l | 


4200 4600 5000 5400 5800 6200 


Fic. 1.—Mean values of A, as a function of wave length 





LOG x, i | | j 
2.00}— k 


0.50}— 


0.00}— 





J ! | | 
-2.00 -1.50 -1.00 -0.50 0.00 = 0.50 








LOG x; 


Fic. 2.—Relation between mean values of log X, and log Xo’ for each Rowland 
£ 


<- 


in tensity from —1 to +6. 
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were plotted one against the other. The straight line passing through 
the points in Figure 2 gives the relation 


log X, = 1.62 + 1.12 log XQ’. (13) 


The equivalent-width calibration is somewhat steeper than that de- 
rived from the multiplet strengths. This result may possibly arise 
from systematic errors in the curve of growth, particularly in the 
flat portion, where, for a given vaiue of log W/X, a small displace- 
ment in the curve may produce a disproportionately large change in 
log Xo. 

TABLE 6 


FINAL SMOOTHED CALIBRATION, CORRECTED FOR OPACITY VARIATION 



































Wave LENGTH 
R | l 
3000 | 3400 3800 4200 4600 5000 5400 5800 6200 6600 
| | | 

+8. 4.21*| 3.47 | 43.17) +2.97\ +2.74)/ +2.71/+2.18)+1.32)+1.57/+1.68 
+7 3.84 | 3.18 |4+2.84/+2.62)/+2.44/-+2.40/+1.92|+1.11]/+1.32/+1.43 
+6 3.45 | 2.89 |+2.50\+2.27/+2 14|/+2 o9|+1 65|+0.90/+1.08)+1.19 
+5. 3.08 | 2.60 |4+2.17/+1.92\+1.84|/+1.78)+1.39|+0.70]+0.83/+0.94 
+4. 2.71 | 2.30 [+1.83)+1.57/+1.54)+1.47|+1.12/-+0.49/-+0 50|+0.70 
+3. 2.34 | 2.02 |-+1.49/+1.22)+1.24/+1.16/+0.85)+0. 29/+0.34/+0.45 
+2 1.95 | 1.73 |+1.16 +0.87|+0.94|+0.85}+0 50|+0.08|+0.09 +o0.21 
+1. 1.59 | 1.43 |+0.82)+0.51\+0.64/-+0.54/+0.32|—0.12|—0.16|/—0.03 
O. 1.22 | 1.14 |+0.49|+0.16)+0 34/0 23|-+0 06) —0.33 —0o.40|—0. 28 
—1. 0.85 | 0.85 |+0.15|—0.20}-+0.05|—0.08]—o 20}—0. 53|—-0.65|—0. 52 
—2. 0.47 | 0.56 [—0.18]—0.55|—0. 25] —0. 39] —0.47|/—0.75|—0.90]—0.77 
—3. | 0.10 | 0.27 |—o si Sie go} —o 50/—0 ie 73° 96|—1.14/—1.01 











* The tabular entries are log Xo. 


As the final step in the calibration, the values of log X%’ in Table 4 
were corrected by the addition of appropriate values of —A, derived 
from equation (12). We have no means of ascertaining whether the 
relation (12) is valid on the violet side of \ 4000, but we adopt this 
assumption for lack of a better one. We then smoothed the values 
of log X, thus derived by the graphical process of plotting R against 
log X, for each wave-length interval separately and representing 
each set of points by a straight line, which was computed by least 
squares. The final calibration is contained in Table 6 and is shown 
graphically in Figure 3. 
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Fic. 3.—Smoothed calibration of Rowland’s scale in terms of log X, as a function 
of R and X. 
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Fic. 4.—Comparison of present calibration with that of Russell, Adams, and Moore. 
Curve A represents log X., averaged over all wave lengths, as a function of Rowland 
intensity for a temperature of 4500°; curve B isa similar plot for a temperature of 5740°; 
curve C shows the variation of log A (determined by Russell, Adams, and Moore) 
with R. The ordinate units of B and C are arbitrary. 
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It is interesting to compare the present calibration with the earlier 
one of Russell, Adams, and Moore. If we plot Rowland intensity 
against mean log X,, averaged over all wave lengths, we obtain a 
curve that may be assumed to correspond with the Russell-Adams- 
Moore plot of R against their log A. The curves are shown in Fig- 
ure 4. Curve A represents the calibration adopted in the present 
investigation with an assumed temperature of 4500°; curve B shows 
a similar calibration for a temperature of 5740°, uncorrected, how- 
ever, for opacity variation, while C is the curve of Russell, Adams, 
and Moore. The ordinate units are arbitrary, since we are concerned 
here only with the slopes of the various curves. Curves A and B are 
considerably less steep than curve C. Also the small difference in 
slope between A and B shows that the type of calibration considered 
here is relatively insensitive to the choice of temperature. 

HARVARD OBSERVATORY 
October 1939 








SOLAR ERUPTIONS 
R. G. GIOVANELLI 


ABSTRACT 


Observations of 268 solar eruptions made with a spectrohelioscope during the period 
June, 1937, to December, 1938, have been analyzed to show the relative distributions 
of intensities and areas. A simple photometer for the measurement of eruption inten- 
sity is described. 

The width of the Ha line during eruptions is constant for small intensities but in- 
creases rapidly when the intensity becomes equal to, and exceeds that of, the continuous 
spectrum. There appears to be no appreciable vertical velocity of the eruption itself. 
Approximately 20 per cent of the eruptions eject eruptive prominences, and, assuming 
that this is caused by La radiation pressure, the width of La is found to be approxi- 
mately six times that of Ha when measured in velocity units. 

The position of eruptions in the chromosphere is discussed, and it is concluded that, 
in general, eruptions are confined to the levels in which the normal bright hydrogen 


occurs. 
OBSERVATIONS 

The observations used for the following discussion have been 
made by the author with the spectrohelioscope of the Common- 
wealth Solar Observatory. Surveys of the sun’s disk have been made 
throughout days of sunshine at intervals of approximately half an 
hour, except during eruptions and periods of considerable solar ac- 
tivity, when the observations were continuous. As far as possible, 
the program of observations following the discovery of an eruption 
consisted of (1) asketch of related spot group and eruption, (2) meas- 
urement of area, (3) estimate of intensity, (4) measurement of width 
of the eruptive Ha line, and (5) examination of the eruptive area 
and its surroundings for velocities of bright or dark hydrogen. These 
were repeated as often as was felt desirable during the eruption. 

The area measurements are affected by the following: (a) poor 
seeing conditions, (6) the narrowness of eruptions, (c) the diffuse 
edges of many eruptions, (d) slight changes of focus owing to tem- 
perature, and (e) the impossibility, under existing arrangements, of 
bringing the measuring graticule into contact with the image. It 
seems that an error of at least 20 per cent may be expected in the 
area estimates. 

The intensity or brightness of an eruption was estimated by eye 
on an arbitrary scale. As a large number of eruptions occurred dur- 
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ing June and July, 1937, it was not long before a fairly consistent 
set of intensities, on a scale from 1 to 5, was being used. Moreover, 
it was soon found possible to subdivide each division into three, e.g., 
2; 2.3; 2.7; 3. Actually the maximum intensity recorded has been 
4.3. Since the intensity was measured by eye on an arbitrary scale, 
it is evident that, over a period of 18 months, without any standard 
of intensity, there were bound to be slight variations, possibly sys- 
tematic, in intensity estimates. Any one observation may have been 
in error by 0.3, or possibly 0.6, units; but it is not considered prob- 
able that an estimate should have been in error by one whole unit. 
The intensity is disturbed by bad seeing or by a lens slightly out of 
focus, as mentioned above, and also by the presence of cirrus clouds 
or sky haze, which reduces the intensity of the image and blurs it by 
scattering. 
TIME OF MAXIMUM INTENSITY 

For twenty-four eruptions the times of commencement and maxi- 
mum intensity were observed. The following summary of these may 
be made. (a) The mean time for an eruption to reach its maximum 
brightness was 7 minutes. (6) For half the eruptions the intervals 
between commencement and maximum were 4 minutes or less. (c) 
The mean interval for eruptions of importance 1 was 4 minutes, as 
against 11 minutes for eruptions of importance 2. It must be remem- 
bered that the mean duration of eruptions of importance 1 is only 
21 minutes, as compared with a mean life of 38 minutes for eruptions 
of importance 2." (d) There was a tendency for the interval to in- 
crease with increased brightness, i.e., the brighter the eruption the 
longer was the time required to reach that brightness. 


DISTRIBUTION OF INTENSITIES AND AREAS 


It is of interest to find how intensities of eruptions vary with 
areas. M. Waldmeier has already stated’ that intensities are more 
or less proportional to areas, but he has not supplied figures to sup- 
port this. In Table 1 the observed intensities are given under ranges 
of the area. It is clear that the intensity increases with the area. It 
is seen that Waldmeier’s statement receives support from these ob- 


t Waldmeier, Zs. f. Ap., 16, 276, 1938. 
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servations, although the variations for individual cases are consid- 
erable. 

In the reports published in the Bulletin for Character Figures of 
Solar Phenomena the importances of eruptions are given on a scale 
of 1 to 3. It was recommended? that, in order to preserve uniformity, 
the “importance” should be determined by its area alone, and for 
these purposes a small eruption was to be called of importance 1, a 


TABLE 1 


INTENSITIES AND AREAS OF ERUPTIONS 





























| | l | 5 
| | | | | | 
Area 25 75 | 125 | 175 | 225 275 325 | 375 | 450 750 All 
a — — meccet cies _ —— | a —_ | — | = - 
Occasions..... | 19 65 30 | & | 14 16 6 | re} 6 | 5 189 
Mean intensity; 1.3 | 3.7 | 2.3 | 2.3 | 2.5 | 2.3] @:4.) 2.7 2.6 | 3.9 | 2.1 
| | | 
| | | | 














Nore.—Areas are given in units of one-millionth of the sun’s disk. 


TABLE 2 


IMPORTANCE AND AREAS OF ERUPTIONS 


! 





| 
een ee 
Number of Waldmeier’s 


| 
Importance Eruptions | Mean Area | 
7 | | Mean Areas 

Used } 

Ae : - eee | 
Ei 78 | 106 | 120 
2. 66 326 380 
cee 4 | 800 | 1020 





medium-sized eruption of importance 2, while very large eruptions 
were to be of importance 3. With the exception of the suggestion that 
eruptions of importance 3 should have an area of 750-millionths of 
the sun’s disk, the estimate of importance was thus left largely to the 
observer. Hence, it is of interest to determine the mean areas for 
eruptions of given importance. These are shown in Table 2, to- 
gether with the only other similar table previously published.’ The 
mean areas given by Waldmeier are about 1.15 times those obtained 
at Mount Stromlo. 

It is desirable to determine the numbers of eruptions occurring 
with different intensities. In Table 3 the eruptions have been sepa- 


2 Bulletin for Character Figures of Solar Phenomena, No. 31, 1935. 
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rated into groups of the same importance, and the number with the 
given observed intensity has been listed. However, as the maximum 
intensity may have occurred some time before the discovery of the 
eruption, this is not a true distribution of intensities. To allow for 
this, suppose provisionally that the intensity of an eruption increases 
linearly to a maximum during a time a and then decreases linearly to 
zero during a time 8. Let i and J be the maximum and the observed 


TABLE 3 


DISTRIBUTION OF INTENSITIES 















































NUMBERS Sicetain 
’ ALCULATED 
NUMBERS OBSERVED OBSERVED MEAN 
. DISTRIBUTION 
(SMOOTHED) 
INTENSITY = 
(/) 
Importance Importance Importance Importance 
Se ete ae 
| 
I 2 3 I 2 I 2 I 2 3 
| | = tele tina ps 
£3. |. I 2 ea I king a Fee rs 2 
4.0 I 2 I I 2.5 4 5 2.5 4 I 
4.9 5 4 I 3 5 Io IO 6.5 7.3 I 
4.2 2 4 5 5.5 13 17 9 rs ae eee 
3.0 9 18 7 10.5 15 16 II a) eee 
2.4 9 3 IO | 10.5 19 II 4.5% 1 JOsG bee dn 
2.3 8 9 | 15 9 27 4 21 | 
2.0 25 12 | 20 7 Pi Bi os a 24.5 IE) A eee 
Se 24 Se fh 22 | 5 re. fei dieok eee st | Seeger 
1.3 16 44 Wf 3 10.5 SD) ee 
1.0 23 4 18 i Pea tavcke esos 9 GSES ncine s 
OF. 10 2 | 10 | | nee e. Poerteretomans 
| | 











| 
| 
| 
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intensities, and f(z) and F(J) the true and the observed numbers of 
eruptions with these intensities. Then it may be shown that 





ee yaw 


where 67 is the smallest measured increment of intensity, i.e., 6z = 3. 
This is the number given in the column ‘Number Observed 
(Smoothed)”’ in Table 3. 

Now, for eruptions of importance 1 mean values are 


a = 4min.; 8 = 17 min. 
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For eruptions of importance 2 
a = II min.; 6 = 27 min. 


The evaluations of f(z) can thus be made, starting, for eruptions of 
importance 1, from J = 4,2 = 0, giving f(4). The value of f(4) may 
then be substituted in the formula for f(3.7), and the process con- 
tinued until f(z) falls to zero. The values obtained are designated 
“Calculated Distribution” in Table 3 and have been reduced so that 
the total calculated number and the observed number are equal. The 
foregoing procedure, however, probably produces an overcorrection. 
If the intensity-time plot had been rectangular, the observed distribu- 
tion would have been the same as the true one. Thus, the actual dis- 
tribution probably lies between the two, and so the mean of these two 
values is adopted and tabulated on the right of Table 3. The note- 
worthy points are: (1) There is a wide range of scatter in intensity 
within any one importance group. (2) The maximum number of oc- 
currences takes place as follows: 


Importance Intensity 
Bene bs Ramee ese ie 2.0 
7 ta eee was 336 
: oe oe 


This shows a progressive increase with importance. 

It has recently been felt desirable to calibrate the scale of inten- 
sities used in the foregoing observations and also to provide a more 
quantitative measurement of the intensity of an eruption. One type 
of photometer for the measurement of intensity with a spectro- 
helioscope has already been described by Newton and Woolley.’ In 
this the intensity of a flocculus is measured by comparing it with a 
photometer field placed immediately above the flocculus. As erup- 
tions frequently have the appearance of bright points or lines on a 
darker background, it seems more satisfactory to compare them 
with a similar photometer image, as described below. 

A isa motor head lamp (Fig. 1), run, with a transformer, from the 
A.C. mains. The light passes through a red filter, B, chosen to match 


3M.N., 96, 5, 1935. 
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Ha, and then through an optical wedge, C. The light is then focused 
by a lens, D; and after passing through a diaphragm, E£, to reduce 
scattered light, the filament is focused onto a slit, F. This cuts off 
all but a very small part of the filament of the lamp. A thin piece of 
unsilvered glass, G, at 45° to the beam reflects a part of the light into 
the eyepiece, H, which is equidistant optically from F and the second 
slit of the spectrohelioscope. Most of the light coming from the sec- 
ond slit passes through G, and the filament image appears super- 
imposed on the solar field. This creates a fairly good “artificial 





Fic. 1.—Schematic diagram of eruption photometer 


eruption” the intensity of which may be varied by using the wedge. 
Newton and Woolley’s method: for determining the scattering from 
one wave length to another is used. Suppose the grating removes the 
fraction a of the light of any particular wave length and scatters it 
in directions proper to other wave lengths; then a may be readily de- 
termined from the known central intensity of Ha (12 per cent), as 
compared with the value measured with the eruption photometer. 
The increase of intensity due to the artificial eruption is measured 
in terms of the continuous spectrum; to obtain the true intensity of 
the flocculus this has to be increased by a factor of 1/(1 — a), since a 
is the fraction of the excess of floccular intensity over Ha which has 
been removed by scattering, and to this has to be added the initial 
corrected intensity of the Ha background, i.e., 12 per cent. 
The optical wedge used obeyed the law 


I = [, 0.10~% 54 , 


where x was measured in millimeters, over about 5 cm and was found 
satisfactory in the range of brightness. The whole photometer is 
extremely simple in operation and, as it reproduces eruption condi- 
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tions fairly satisfactorily, may be recommended with confidence for 
eruption photometry. 

For the calibration of intensities the intensity of an eruption was 
first estimated by eye and then measured with the photometer, the 
intensity of the continuous spectrum background 2.6 A from the 
center of Ha being determined immediately afterward. Not many 
observations have yet been made, but it has been thought worth 
while to include them to give a quantitative idea of the foregoing 
scale of intensities. Mean values are given in Table 4. 


TABLE 4 


CALIBRATION OF INTENSITY SCALE 








eee pe ere eee ee on, 4 21 3 
Percentage of continuous spectrum....... 36 70) II 





WIDTH OF THE ERUPTIVE Ha LINE 

On many occasions during eruptions the width of the Ha line has 
been measured. To do this the line-shifter was rotated toward the 
violet until the eruption just disappeared, and the position was 
noted. The same procedure was followed on the red side of the line, 
and the position of the center of Ha was also obtained. The scale 
being calibrated in units of kilometers per second, the width was also 
given in terms of that unit, although whether or not the entire width 
is due to Doppler effect need not concern us here. The settings of 
the line-shifter were usually reproducible to within 3 km/sec. In 
nearly all cases the intensity of the eruption was estimated at the 
same time. The widths of Ha are grouped under intensities in Ta- 
ble 5. It will be noticed that the width, which is more or less con 
stant for values of J up to J = 2.3, begins to increase for J = 2.7 
(which corresponds to about 105 per cent of the continuous spec- 
trum) and increases rapidly thereafter. The increase thus begins 
when the central intensity becomes approximately equal to the 
neighboring continuous spectrum. 

VELOCITIES OF ERUPTIONS 


It is of importance to decide whether an eruption is stationary in 
the chromosphere or whether it rises, as do eruptive prominences. 
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During all the author’s observations no eruption has been ob- 
served on the disk in which there was a definite displacement by 
velocity of the position of maximum intensity. Moreover, by the 
courtesy of Dr. C. W. Allen, the author has been permitted to ex- 
amine several spectra and microphotometer tracings of eruptions 
made by him at the observatory on which also, although there is a 
dark reversal in the center of the bright line, it is quite clear that 
during these eruptions there was no appreciable displacement of the 
eruptive line. Actually, the dark reversal produces a flatter line and 


TABLE 5 


WIDTH OF ERUPTIVE Ha LINE 





























= : | 
Intensity | I r.3 | 1-7 2 2.3 2.7 3 3.3 3-7 4 4-3 
=e ass ea seein nm can! aM a ala a a 
Number of ob-| 
servations... 9 10 | 23 21 16 14 20 II 8 4 I 
- 
Mean width | | 
(km/sec).....| 82 74 | 86 84] 85 | QI | 105 | 123 | 124 1g2 
| | } | 
| | | } | 








makes accurate visual measurements of any very small displace- 
ments difficult. Further evidence on the stationary state of erup- 
tions is obtained from the width observations previously mentioned. 
If we assume that the emission line is symmetrical—in the wings, at 
any rate—then the mean of the blue and red edges of Ha gives a 
measure of the velocity. The distribution of velocities as determined 
in this way is given in Table 6, where the probable error of a deter- 
mination is about 3 km/sec. The results vary considerably but 
show a recessive mean of 3.4 km/sec. If this represents a true veloc- 
ity, it indicates that any velocity is, on the average, inward. 

On one or two occasions the brightest part of the eruption was not 
observed quite at the center of Ha; on these occasions, however, it 
was clearly seen that a dark marking, due to an overlying promi- 
nence, was superimposed on the eruption; and because of a vertical 
velocity, its absorption line was removed to one of the sides of Ha, 
thus distorting the eruptive line. Perhaps where displacements of the 
Ha emission line have been previously reported, phenomena of this 
type may have occurred. 
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OBSERVATIONS OF OTHER LINES 






The following lines were, at one time or other, observed to be 


affected during eruptions: 


a a Ha HB 
NES Seen Sante es D, D, 
Magnesium........... b, b2 

Iron (neutral)......... 5269.7 

Iron (ionized)......... A5018.6 b, 

aC re! 3 


Hy 


b, 


In all elements other than hydrogen the effect during eruptions is 
very weak, and the observations are most incomplete. Even so, a 


TABLE 


6 


VELOCITIES OF ERUPTIONS 




















Velocity Number of “ 
¥ e Total 
(Km/Sec) Occurrences 
II I 
Q -II ° 
Approach... { 64- 83 I 35 approach 
Chae 5 
\ 13- 34 12 
—-I-1 30 II zero 
' 19- 33 26 
4-6 24 
6}3- 8} 15 
‘ Q -II ‘ 
Recessive. . { 9, 1 3 g2 recessive 
113-133 2 
14 -10 4 
164-18} 3 
| \ 19 3 
! 
Mean of approach velocities........ — 2.1 km/sec 
Mean of recessive velocities. . +5.9 km/sec 
Bi OG ty VOCE oc aio. i5 ae seo cuales +3.4 km/sec 


few generalizations may be made. (1) In all elements except helium 
an increase in radiation was observed. With the D, helium line, not 
usually visible in the solar spectrum, the effect produced was an ab- 
sorption. (2) The effects in the foregoing lines can be observed only 
when the Ha intensity is high. (3) The position and shape of the 
eruption in other lines is that of the brightest parts of Ha. (4) The 
eruptions can be observed only in the very center of the line. Since 
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these lines are all much narrower than Ha, any slight velocity would 
have displaced the emission into the wings or even into the continu- 
ous spectrum. From (4) it is clear that, in these cases at least, the 
elements concerned had no appreciable vertical velocity. 


ERUPTIONS AND ERUPTIVE PROMINENCES 

In view of the foregoing conclusions it is of interest to review opin- 
ions on the nature of eruptive prominences. Of early workers, Young 
has said: 


In a few instances the gaseous eruptions in the neighbourhood of a spot are 
so powerful and brilliant that with a spectroscope their forms can be made out 
on the background of the solar surface in the same way that the prominences are 
seen at the edge of the sun. In fact, there is probably no difference at all in the 
phenomena, except that only prominences of most unusual brightness can thus 
be detected on the solar surface.4 


Abetti expressed the opinion’ that “the bright flocculi are, as a rule, 
eruptive prominences,”’ while McCrea® has recently come to a sim- 
ilar conclusion. Fox,’ on the other hand, stated that eruptions were 
the bases of eruptive prominences. Pettit remarked* that, in gen- 
eral, no markings appear on disk plates (presumably spectrohelio- 
grams) in the region of eruptive prominences and, with McMath, 
further supported this view? in objecting to Kiepenheuer’s theory of 
eruptive prominences, pointing out that eruptive prominences may 
be widely distributed over the disk, while eruptions are limited to 
spot zones. 

Since eruptions have no vertical velocities, the remarks of Young, 
Abetti, and McCrea seem, in general, to be without support. It will 
be shown below that there might be exceptional occasions, however, 
in which their views are valid. With the other observers the differ- 
ences in opinion lie in confusing two distinct types of eruptive 
prominence. (a) The first type commences as a quiet prominence and 
for some, as yet unknown, reason acquires a velocity outward and 
rapidly disintegrates. This is the type most frequently examined by 
Pettit. (b) The second type is that described by Fox, i.e., a type as- 

4 The Sun, p. 135, rev. ed., 1904. 


5 Handb.d. Ap., 4, 121. 
6 Observatory, 62, 16, 1939. 


7 Ap. J., 28, 253, 1908. 
8 Pub. Yerkes Obs., 3, 205. 
9 Ap. J., 88, 244, 1938. 
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sociated with eruptions. It is possible that some of these correspond 
to those listed by Pettit as ‘‘surges.” 

When observed on the limb, an eruptive prominence shoots out 
with high velocity and in the course of an hour usually disintegrates, 
but at times parts of the prominence may be seen turning over and 
falling back into the sun. In twenty cases since the beginning of 
March, 1938, when this matter was first given attention, eruptions 
have been seen on the disk, and strong absorption markings with 
high approach velocities have been found in the neighborhood. These 
markings are of a quite different type from those described by Hale,’ 
which usually show approach velocities at one end and recessive 
velocities at the other and take the form of neat curved lines. The 
dark markings referred to above, however, occupy a wide irregular 
area. Usually they are first seen devoid of accompanying recessive 
dark markings; but the latter, equally irregular in shape, soon ap- 
pear, while the approach markings decrease in intensity and disap- 
pear after an interval of about half an hour. It is quite apparent that 
these masses of dark hydrogen are really eruptive prominences, 
ejected from the eruption proper and strongly absorbing the wave 
lengths corresponding to the Doppler velocities. This view is com- 
pletely verified by the observation of five of these eruptions on the 
disk but near the limb, in which the dark hydrogen showing veloc- 
ities traveled to and passed beyond the limb, continuing outward to 
appear as eruptive prominences. Additional evidence, if necessary, 
was obtained in the case of the eruption of December 19, 1938 
(to be discussed later), which was situated right on the limb of the 
sun. Although not in vertical motion itself, a prominence was ob- 
served to be expelled from a neighboring area and to reach a height of 
65,000 km before disintegrating. In ten further cases approaching 
dark flocculi, obviously eruptive prominences in projection, were ob- 
served without any eruption but near active spots. In one of these 
it was subsequently discovered that a radio fadeout had occurred 25 
minutes previously, indicating that an eruption had probably taken 
place around this center. It seems that in these cases the eruption 
had died down before the observations were made, leaving only the 
eruptive prominence projected on the disk. 


10 Ad. J ., 71, 55; 1930- 
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It is quite conclusive that this second type of eruptive prominence 
does exist and is expelled from eruptions. During the period March-— 
December, 1938, 112 eruptions were observed, of which 21, or 19 per 
cent, gave rise to eruptive prominences. What, then, is the cause of 


TABLE 7 
ERUPTIONS PRODUCING ERUPTIVE PROMINENCES IN PROJECTION 



































VELOCITY OF 
PROMINENCE 
DATE bapa arom AREA seiamell ere REMARKS 
(E.A.S.T.) TANCE a il SITY 
Mean Mant 
mum 
1938 
Mar. 17 1254-1256 I rel 48 2,7 al Ee ee rer 
Apr. 11 1538-1546 I een ee 97 | Prominence ex- 
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June 1 1105-1147 I 3a beers 87 Same 
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Dec. 7 | 1247-1250 | I | 3-3 |------ 330 Same 
| | | 








this special phenomenon accompanying some eruptions? It might 
have been expected that the chance of producing eruptive promi- 
nences would have been greater for the eruptions of greater inten- 
sity. In Table 7 are given the eruptions, together with area, impor- 
tance, intensity, and velocity, of the accompanying prominence. The 
mean area is 175, compared with a mean for all eruptions of 188, 
showing that the chance of producing an eruptive prominence is in- 
dependent of the importance, or area, of the eruption. Again, the 
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mean intensity is 2.0, compared with a mean observed intensity of 
2.1 for all eruptions. None of these eruptions was observed at its 
beginning, but that does not seem sufficient to account for this re- 
sult: the chance of producing a prominence is independent of area 
and intensity. Remembering that the other type of eruptive promi- 
nence arises from an existing prominence, it seems probable that the 
presence of a low-lying prominence is essential before the second type 
of eruptive prominence can be produced, and that when the condi- 
tions of position of prominence and eruption are suitable, the erup- 
tion will give rise to an eruptive prominence. It is rare for an erup- 
tive prominence to arise from the eruption itself, but there is usually 
a gap of about 10,000 km between the dark markings and the erup- 
tions. 

The prominence velocities given in Table 7 represent the resultant 
of Doppler velocity measurements and movement at right angles to 
the line of sight. Let us see if the velocity of the prominence is af- 
fected by area or intensity of the eruption. The mean of the maxi- 
mum velocities of prominences associated with eruptions of impor- 
tance 1 is 187 km/sec, while with eruptions of importance 2 it is 
153 km/sec, a result showing no appreciable effect of area. However, 
when the maximum velocities are compared with //a intensity, there 
is a definite tendency toward an increase in velocity with intensity. 
If the velocity of 320 km/sec, occurring on November 24, 1938, beomit- 
ted on the assumption that the eruption had been of greater intensity 
before it was first observed, the velocities fall into two groups: (a) 
small velocities (average, 114 km/sec) for intensities 2.3 and less; 
(b) high velocities (average, 271 km/sec) for intensities 2.7 and more. 
This is exactly the same division of intensities as has been found to 
influence the width of Ha. It is thus probable that these prominences 
are expelled by the mechanism of Kiepenheuer’s theory," i.e., the 
light-pressure of Lyman radiation from the eruption, for with greater 
intensity, the width of Ha, and presumably La, would be greater, 
and consequently the velocity of the prominence would be high. 

If it be accepted that these prominences are ejected by radiation 
pressure of La, the observations may be used to estimate the relative 
widths of Ha and La. However, as in many of these cases the maxi- 


11Zs. f. Ap., 15) §3) 1937: 
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mum intensities were not observed and in most of the Ha widths 
were not determined, this estimate can only be made indirectly. 

The distribution of intensities among eruptions has been given in 
Table 3; and, by using the mean widths of Ha given in Table 5, the 
relative numbers of eruptions expected to have maximum Ha widths 
in given ranges were obtained. The distribution of prominence 
velocities was smoothed for comparison purposes. Now, if to a given 
Ha width there always corresponds a definite La width, if we neglect 
the inverse-square effect, and if, in the foregoing cases, there has been 
no other mechanism acting to increase the velocities of the promi- 
nences, then, corresponding to a given Ha width, there will be a defi- 
nite prominence velocity. By a direct comparison between the two 
foregoing distributions this relationship has been determined. It 
may be expressed by 


Velocity = 6.2 (} Ha width — 25) , 


when the Ha width is expressed in velocity units. 

The width of La, then, in velocity units is approximately six 
times that of Ha. It may be pointed out that the theory of radio 
fadeouts advanced by Martyn, Munro, Higgs, and Williams” re- 
quired a width of La of approximately this amount. The constant in 
the preceding expression would be interpreted as representing the 
distance from the edge of the La line at which radiation pressure is 
just sufficient to balance gravity. 


POSITION OF ERUPTIONS IN THE CHROMOSPHERE 


It has previously been pointed out'’ that, where the base of an 
eruption is obscured by the limb, radio fadeouts have not been ob- 
served, while fadeouts can and do occur in association with eruptive 
prominences whose bases are on the visible side of the limb. This is 
taken to indicate that the eruptions causing these prominences are 
at no great height above the chromosphere. 

Again, on several occasions eruptions have been observed pro- 
jecting over both umbra and penumbra of sunspots. There has been 
no detectable discontinuity or change in the appearance of the erup- 


12 Nature, 140, 603, 1937. '3 Giovanelli, Ap. J., 88, 204, 1938. 











348 R. G. GIOVANELLI 


tions at these places, and hence the eruption lies at a higher level 
than umbra or penumbra. 

To confirm these statements, a watch was kept for an event on the 
edge of the sun which could unhesitatingly be classed as an eruption 
proper on the limb. Such an eruption was observed on December 
19, 1938. In height, the eruption was 10,000 km, half of which pro- 
jected beyond the edge of the hydrogen chromosphere; this is con- 
siderably greater than the usual width of 3,000 km. Now, the height 
of the chromosphere is about 8,000 km, so that the center of the erup- 
tion was about 8,000 km above the reversing layer. The intensity of 
the eruption, which seemed uniform throughout its height, bright- 
ened considerably from 11"05™, when first observed, to a maximum 
value of 3 at 11508™. By 1151g™ the intensity had decreased to 1.7. 
The eruption itself showed no tendency whatever to rise, but at 
11'97™ there was observed the ejection of an eruptive prominence 
fainter than the normal Ha disk, which reached a height of 65,000 
km. The base of this prominence did not quite touch the eruption, a 
gap of about 10,000 km separating them. The whole sequence of 
events was so similar to that observed during an eruption on the disk 
that it was felt that no mistake in identity could have been made. 

Since eruptions generally appear as brightenings in already bright 
hydrogen flocculi, and as the latter usually occur in or just above the 
normal hydrogen chromosphere, we may derive the following con- 
clusions about the general nature of eruptions. (1) Eruptions take 
place in bright hydrogen flocculi. As they show no vertical velocity, 
they do not leave their original position during their existence. (2) 
Approximately one-fifth of the eruptions give rise to eruptive prom- 
inences whose intensities are less than that of the center of Ha. 
These do not usually touch the eruption itself but are separated by a 
gap of 10,000 km or more. If these be ejected by radiation pressure 
of La, it appears that the width of La must be about six times that of 
Ha, when measured in velocity units. (3) The brightness in Ha is of 
the same order as the continuous spectrum in the neighborhood of 
Ha. The width of the eruptive Ha line, which is more or less con- 
stant for lower intensities, increases rapidly when the intensity be- 
comes equal to and exceeds that of the continuous spectrum. 

This discussion seems to give fairly definite conclusions as to the 
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general conditions of velocities and positions of the eruptions in the 
chromosphere. However, it must be realized that exceptions are 
possible. To illustrate this, on March 14, 1939, at 10"15™ E.A.S.T., 
there was observed to grow, most probably from an active group 
which had then passed 3° behind the limb, a brilliant prominence 
whose maximum intensity reached go per cent of the continuous 
spectrum. It reached a maximum height of 84,000 km at 10%33™ and 
attained a maximum area equal to 4oo-millionths of the sun’s disk. 
By 10"42™ it had shrunk and completely disappeared. In all respects 
it seemed similar to the great surges described by McMath and Pet- 
tit. This prominence must surely have been classed as an eruption 
had it been situated on the sun’s disk, for the intensity was equal to 
that of the average eruption. According to a recently developed 
formula,"’ an eruption of such intensity and area would not have been 
expected to produce a complete radio fadeout; yet it would probably 
have caused a weakening in £-echoes; and this is precisely what hap- 
pened. It must be concluded, then, that this eruptive prominence, 
or surge, was indeed an eruption. Unless the apparent rise of the 
eruption was due to an upward transfer of luminescence—a phenom- 
enon whose existence is at present doubtful—this represents an ex- 
ception to the general rule that eruptions have little or no vertical 
velocity, for the mean velocity recorded was 78 km/sec. From the 
evidence previously presented, it must be deduced that such erup- 
tions are rare. 
COMMONWEALTH SOLAR OBSERVATORY 


CANBERRA, AUSTRALIA 
July 1939 


14 Ap. J., 85, 296, 1937. 's Giovanelli and Higgs, Terr. Mag., 44, 181, 1930. 











SOME CORRECTIONS TO DREYER’S CATALOGUES 
OF NEBULAE AND CLUSTERS* 





DOROTHY CARLSON 


ABSTRACT 


Corrections to Dreyer’s catalogues compiled from Mount Wilson plates and other 
photographic sources; 330 objects are to be dropped from the New General Catalogue 
(Table 1a) and 250 from the two Index Catalogues (Table 1b). Positions of 49 objects 
are corrected (Table 3), and the existence of 23 doubtful ones is confirmed (Table 4). 





The Mount Wilson collection of photographs made with the 1oo- 
inch and 60-inch reflectors includes about two thousand extra- 
galactic nebulae occurring in Dreyer’s New General Catalogue and 
about seven hundred in the /ndex Catalogues, in addition to a con- | 
siderable number of galactic nebulae and clusters. In the course of 
identifying and classifying these objects, various lists of corrections 
to Dreyer’s catalogues have been consulted, and a new list has been 
compiled from the Mount Wilson plates. The publishing of the 
Mount Wilson list offers an opportunity to bring together the more 
authentic corrections from all sources. 

Experience demonstrates the great reliability of corrections de- 
rived from photographs in comparison with those derived from 
visual observations. Consequently, the corrections included in the | 
present list are, in general, restricted to photographic sources. How- | 
ever, Dreyer’s own corrections, assembled for the most part from 
visual observations, have been freely consulted; and a supplemen- 
tary list (Table 2) has been formulated containing objects which 
Dreyer proposed to remove from the catalogues, although no photo- 
graphic evidence is available for confirmation. The lists of correc- 
tions, although extensive, are far from complete, and additions will 
accumulate for a long time to come. 

The sources for the present compilation, together with the sym- 
bols by which they are represented in the tables, are as follows: 





W Mount Wilson photographs. 
L, Publications of the Lick Observatory, 13, 1918. 


* Contributions from the Mount Wilson Observatory, Carnegie Institution of Wash- | 


ington, No. 626. 
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R Heidelberg Veréffentlichungen, 9, 1926. Reinmuth attempted to describe 
the 4445 objects north of Decl. = —20° in John Herschel’s general cata- 
logue, as they appeared on camera plates or, occasionally, on plates taken 
with the 27-inch reflector. Uncertainties, in general, are restricted to 
small, compact nebulae described as nebulous stars, for which the small 
scale employed was not entirely satisfactory. The effects are illustrated 
in Table 4, which contains nebulae not found by Reinmuth but found on 
Mount Wilson plates in their catalogued positions. 

Hel Helwan Observatory Bulletins, Nos. 9, 15, 21, 22, and 30. Reports of a 
survey of nebulae between Decl. o° and —45°, made with the Helwan 
reflector. 

H Harvard Annals, 88, No. 1, 1930. A catalogue of 2778 nebulae in the 
region of the Virgo Cluster; includes a list of 39 identical nebulae and one 
of 94 missing NGC and IC objects. 

DI Notes following the NGC, Memoirs R.A.S., 49-50, 212, 1888. 

D II Notes and corrections to the NGC (following First Index Catalogue), 
Memoirs R.A.S., §1, 225, 1895. 

D III Notes and corrections to the NGC (following Second Index Catalogue), 
Memoirs R.A.S., 59, 186, 1908. 

DIV Dreyer’s own corrections to the NGC resulting from the revision of Sir 
William Herschel’s three catalogues of nebulae, Memoirs R.A.S., 73, 37, 
1912-13. 


The other symbols in the tables are those used by Dreyer with the 
addition of “nf,’’ meaning ‘‘not found.’”’ Tables 1a and 16 contain 
listed numbers which, on the evidence of photographic plates, should 
be removed from the New General Catalogue and the Index Cata- 
logues, respectively. Table 2 is a compilation, from Dreyer’s notes, 
of missing nebulae which have not been checked by photography. 
Table 3 gives a few corrections to catalogue positions where the 
errors are more than 5’; and Table 4, nebulae not found by Rein- 
muth and other observers, but which have been found on Mount 
Wilson plates in the catalogued positions. 











































































NGC 
18 
3° 
46 


82 
84 
go 
153 
156 
158 
160 


747 
760 
763 
771 
843 
952 
983 
1040 
1059 
II4I 


1210 
1240 


1142/ 


TABLE 1a 





Osyects To BE STRICKEN FROM THE New General Catalogue 


aad W 

**W,R 

ae 

No trace of nebulosity, W, L. 
“VF, pL, pmE, ph, = neb Wolf 
XIII No. 112; neb Wolf XIII No. 
111 np 9o’,” R. However, Wolf 
XIII No. 112 is 1™6 f and 5's of 
Dreyer position of 56 and is pL 
(1/7 in diam) instead of eL as 
Herschel has described it. 

*W 

*W 

*W 

= NGC 151, W, DII 

** W 

nf, W 

= NGC 162, W. In NGC, h39 sp 
should read h32 sp. NGC 162 is 
merely a better position of NGC 
160 and is listed as such by 
Schultz. R mistakes a star for 
162. 

nf, W, DIV, R 

“Nothing shown here, but 4’ s of 
IC 41 there isa neb, pF’, S, E go°,” 
Hel 

** R 

4 st, 10-11.5, in parallelogram, R 
nf, W,R 

*** W, R 

OW UR 

nf, R 

nf, plate not good, R 

= IC 2107, RK 

= iC 1117, R 

eet WR, DEL) 

ee W, D II, R, Hel 

ae R 

wk. aT. Dal 

= NGC 6977, R 

*** R 

*W 

FW -DILR 

nf, Hel 

nf, Hel 

ae W 

nf, Hel 

50 Cas, nebulous?, R 

*** R 

nf, R 

nf, R 

= NGC 1053, R 

**R DI, DII 


**R 


nf, Hel 
RR 


NGC 
1251 
1209 
1276 
1312 
1318 
1340 
1307 
1369 
1378 
1429 
1436 
1442 
1443 
1458 
1488 
I551 
1724 
1746 
1750 


1758) 


4o/ 
1781 
1908 
1909 
1927 
1988 
IQQO 


2045 
2001 
2107 


2189 
2330 
2334 
2301 
2386 
2387 
2390 
2391 
2404 
2471 
2492 
2515 
2542 
2597 
2043 
2053 
2703 
2705 
2707 
2727 


** R 

nf, Hel 

** W 

oP ae 

nf, W, Hel 

= NGC 1344, W, DIII 

= NGC 1371, W 

nf, W, Hel 

** W, Hel 

*W 

nf, Hel 

= NGC 1440, R, DI, DIV 
nf, W 

nf vis, Hel 

** R 

= NGC 1550, 
*e* R 


R, Di, DIV 


vL Cl, R 


of. Rk. DID 
= NGC 1794, R, D III 


nf, R 

nf, R 

df. RDI 
2. Dy 


nf, by Curtis on 2 hr. exp. (L, 13), 
who suggests that Herschel was 
here misled by the radiance about 
the bright star, or that he ob- 
served the similar nebulosity,IC 
434, south of ¢ Ori instead. 
BD+12°884, no nebulosity,W, R 
No Cl, Hel 


“No nebulosity found around 
BD—6°1412,” Hel 
nf, R 


= IC 4c8?, R, DN, DIV 
= IC 464?, R, DIT, DIV 
= NGC 2359, Hel 

*** W 

‘rz 


** R 

No nebulosity, W 

nf, R 

== IC 2390, R 

gp of st, W 

nf, Hel 

nf, Hel 

nf, Hel 

“No nebula in this position which 
is near BD —2°2766,” Hel 


nf, W,R 

















TABLE 


nf, R 

nf, R 

2 neb and *, W 

nf, Hel 

gp of st, W 

nf, W, Hel 

nf, W 

= IC ss6?, R 

WwW, 

** W 

* W 

** R 

nf, R 

nf, R, DI 

BD+sr 1585, R, L 

= NGC 3144, R, DIV 
= NGC 3155, R, DIV 
i 

= NGC 3183, R, DI, DIV 
nf, R, DI 

*R, DIV 

*W 

** W, Hel 

* 


R 
**W DI,DIV 
*W,R 
*W,R 
= NGC 3329, R, DI, DIV 
nf, R, DIV 


nf, W 

= NGC 3528, R, DIV 
**“** R 

= NGC 346s, R, DIV 
= NGC 3571, Hel 
*R 

nf, R 


«= NGC 3626, W, R, DIV 
= NGC 3645, W, R 

Nest of nebulae, W 
* 


xC 3732, W 

iC 3301, W, R, DI, DII 
1C 3786, W 
x( 
i€ 


tou il 


> 3788, W 
> 3794, W, R, DIV 


il 
ee ae Ot ee te, 


AO RRR EREREZELE 


> 3848, R 


. =] 
* | =n 


ll 


= 
a 4 
~~ 


> 3912, R 


es 53 3 
abe 
O 


3C 39077, W 
, DIII 
NG C 4005, W, R, DIV 


ll ell 


= 
AWA RZ 


ll 


GC 4032?, R, DII 


1a—Continued 


NGC 
4955 


4099 
4107 
4113 
4119 
4139 
4140 
4153 
4154 
4167 
4208 
4209 
4223 
4228 
4205 
4297 
4315 
4317 
4322 
4327 
4338 
4345 
4354 
4357 
4304 
43006 
4397 
4368 
4397 
4398 
44°97 
4409 
4420 
4427 
4437 
4443 
4447 
4471 
4595 
4508 
4512 
4549 
4554 
4500 
4577 
4582 
4004 
4010 
4024 
4029 
4037 
4004 
4667 
4700 
4729 
473° 
4752 
4766 
4797 


353 


“pF-cF, cS, iR, glbM; in Dreyer’s 
place no pB neb found,” R 

nf, R, DIV 

*R,H, DIII 


nf, R, D III 

nf, R, H, D Ill 

nf, R 

nf, R 

nf, R 

- 

= NGC 4163?, 

= NGC 4212, W. R, H, DIV 
nf, W "Rk 

= NGC 4241, R, H, DIV 
= NGC 4214, W, R, DI, L 
= NGC 4263, W, DIII 

nf, R, H 

nf, W, H 

— 


NGC 4323, W, H 


*1 Bil 


IC 3247, R 
W : 


351, W, H 
381, R, DIV 


re hey 


PW, D III, H 
*W, DIT, R, H 
= NGC 4413, W, H, L 
= NGC 4420, W,h, R 


"i, Ree 


= NGC 4417, W, R 
*W. Dill 2..o 

= NGC 4446, W, H 
*W_ ko 

= NGC 4496, W 
nf, H, R 

= NGC 4521?, R 
nf, R 

nf, H, D Il 

nf, R, DIII 

= NGC 4543?, R 
*R 


nf, Hel 

= NGC 4470, W, R, DIV 
= NGC 4636?,R 

** R 


nf, H; = NGC 4647?, R 

= NGC 4665, W, DIV, R 
nf, R, H, DI, DIT, L 
nf, Hel 

nf, Hel 

nf, Hel 

mie 

7er WL 

nf, R, DUI 








NGC 
4804 
4882 
4884 
4910 
4913 
4960 
4962 
4979 
5003 
5008 
50607 
5106 
5111 
5113 
5160 
5200 
5242 
5268 
5310 
5317 
539° 
5404 
5428 
5429 
5432 
5435 
5407 
5527 
5588 
5648 
5699 
5793 
5794 
5819 
5841 


5871 
5877 
6002 
6049 
6064 
6125 
6199 
6354 
6430 
6431 
6466 
6481 
6499 
6550 


6625 
6647 
6660 
6682 
6724 
6728 
6737 










IRREP ON I 


#5 5 
# Fh eh 


NGC 4802, R, Hel 
NGC 4886, W 
NGC 4889, W, R 


,R 

aes 

on 

-R, DIV 
IC 4108, R, DIV 
,R 
Ae °4 
’R 


pes Ey 


nf, 


** 
a 
me 
ee 


=) 
. 


nf, 
nf, 


R 
R 


NGC 5110, R 


GC s1o9, R, DI, DIV 


N 
N 

R, DI 
R 

W 

R 

= NGC 5371, W,R 
*R 


W, Hel 
W, Hel 
W, Hel 
W, Hel 


NGC 5706, R, DIV 
NGC 5709, R, DIV 


R 
NGC 5808, R 
W,R 

W, DIV 


NGC 6052, R, DIV 
NGC 6130, R 


R 
4 st, Hel 


chain of 4 st, R 


nf, 


R 
NGC 6478, W 


*W 


** 


R, L 


NGC 6548, W, D III 


No Cl,R 
No Cl, R 
No Cl, R 


NGC 6661, W, DII, R 


No Cl, R 
No Cl, R 
No Cl, R 
No Cl, R 


TABLE 1la—Continued 


NGC 
6749 
6773 


6795 
6828 


6837 
6839 
6843 
6847 
6873 
6895 
6896 
6901 
6904 
6933 
6947 
6952 
6965 
6966 
7011 
7028 
7943 
7071 
7122 
7150 
7101 
7186 
7202 
7210 
7234 
7255 
7320 
7327 
7333 
7338 
735° 
7403 
7431 
7438 
7439 
7447 
7472 


7565 


I~ 
moe 


NNN Ss 
Dunne 
3 
Wn 


On 
= 
> 


No Cl, R 
No Cl, R 
No Cl, R 
No Cl, R 
No Cl, R 
No Cl,R 
No Cl, R 
nf, R, DI 
** R 


No Cl, R 

nf, R 

= IC s5000?, R 
No Cl, R 

OW, BL 

nf vis. (twice), Hel 
= NGC 6951, W, 
* W, R, Hel 

** W, Hel 

No Cl, R 

nf, > 

ni, R 


= NGC 7067?, R 
*R 


D Ill 


nf, R 

** R 

7Fst,R 

nf, Hel 

nf, R, DI 

No Cl, R 

nf vis. (twice), Hel 
sea,’ Me F 


*R, DI, D II 
nf, R, DII 
4st, R, DII 
72 


3 F st inline, R 
nf, R, DI 


**R, DI, DII 
nf, R 






































| | 
I 
39 
47 
89 
117 
IQ! 
280 
281 
324 
339 
394 
487 
532 
587 
629 
684 
689 
795 
788 
793 
805 
811 
820 
839 
847 
958 
973 
974 
1026 
1098 
1126 
1148 
E72 
1175 
1179 
1184 
I1gO 
1282 
1300 
1307 
1312 
1310 


1325 
1320 
13603 
I44I 
1450 
1457 
I511 
1538 
1539 
1547 
1556 
1658 
1686 


TABLE 10 


Osjects To BE STRICKEN FROM THE Index Catalogues 


**K W 
NGC 178, Hel 

Hel 

"NGC 446, R 
= NGC 558, W 

= NGC 794, R 
gp of st, W 

= NGC 1177, W,R 
= NGC 1331, W, Hel 
* W, Hel 
nf, W 

= NGC 2494, R 
nf, Hel 

BD— 1°2334, W 

= NGC 3312, W 

= NGC 3644, W, R 
= NGC 3661, R 

wK W 
NGC 4405?, W, H 
NGC 4445, W, R, H 
NGC 4611, W, H 
NGC 4663, W 
NGC 4676?, R 
VGC 4851?, R 
NGC 4974, R, D Ill 
nf, W 
nf, W 
nf, W 

*W 

*W 
NGC 5952, R 
NGC 6020?, R 

nf, W 
nf, W 

= - NG C 6054, W 


. Ell 


SEA ES: 


run tu 


0 ll 


= NGC 6892?, R 

= IC 5000; Decl. of 
wrong, of 5000 correct, W 
= NGC 6927, W,R 

= NGC 6928, W 


1316 is 


> 7240, W 


= NGC 7767, R 
> 70, W, R 


= NGC 444, R 
= NGC 4099, R 


wn 


IC 

1710 
1712 
1751 
1765 
1766 
1851 
1871 
1983 
2040 
2041 
2048 
2077 
2091 
2092 
2107 
2133 
2154 
2173 
2232 
2262! 
2266 | 
2270} 
2272| 
2273 | 
2274 | 
2275 
2276) 


2277} 


2278 | 
2279) 
2280 | 
2281 | 
2284 
2289 | 
2290 
2291 
2408 
2410 
2413 
2424 
2425 
2460 
2500 
2509 
2528 
2571 
2579 
2005 
2009 
2610 
2011 
2658 
2659 
2662 
2063 
2064 


NGC 575, R 
NGC 584, W, R 
NGC 741, R 
NGC 783, R 
NGC 785, W, R 
nf, 15 so™ exp., L 
nf, W 

nf, W 

nf, Hel 

nf, Hel 

nf, Hel 

= NGC 1608, R 
gp of F st, W 

** W 

NGC 1707, R 


W 
NGC 2291, W 
NGC 2543, W, R 


iu wt ue 


nf, 
** 


nf, L. Curtis remarks (L, 13): ‘‘I 
am unable to identify any of these 
objects. IC 2274 appears to be 3 
small stars almost in contact. 
Most of the others appear to be 
simply very faint stars. The re- 
gion was carefully checked.” 


*W 
= NGC 2667, W, R 
*W 


= NGC 2704, R, DIV 
nf, Hel 

nf, W 

nf, W 

*W 

+ W 

= NGC 3223, W, Hel 
= NGC 3251, R, DIII 
nf, W 

? = NGC 3404, Hel 
nf, W 


*W 
bbl. 
vil 
7 
nf, W 
*W 








IC 

2726 
2728 
2731 
2865 
2866 
2868 
2970 
3009 
3011 
3018 


3026 
3927 


3°35 
3042 
3°45 
3048 
395° 
3951 
3050 


3959 
3004 
3070 
3071 
3072 
3070 
3083 
3085 
3086 
3087 
3088 
3090 
3098 
3102 
3103 
3113 
3114 
3117 
3123 
3124 
3129 
3130 
3132 
3133 
3139 
315° 
3158 
3160 


3101\ 
3162) 


3103 
3181 
3182 
3183 
3190 
3191 


TABLE 1 


*W 

*W 

*W 

*W 

*W 

= IC 698, W 

ok Ok Ok W 

nf, H 

= NGC 4124, W, H 

nf, may be neb 30’ s of Frost’s 
published position for IC 3018, H 
*W 

nf, a defect on the plate where it 
was originally found by Frost, H 
= NGC 416s, W, H 

= NGC 4178, W, H 

nf, H 

7 OW, 

= NGC 4189, W, H 

= NGC 4193, W, H 

nf, may be neb 1° of Frost’s pub- 
lished position for IC 3056, H 
nf, W 

= NGC 4206, W, H 

"AW, 

nf, H 

nf, H 

nf, H 

nf, H 

*W, H 


= NGC 4241, H 
*W,H 
NGC 4246, W, H 


’ 


H 
H 


**# #5 5 
—— 
=e 


ey 
vot be 


rw", 


wie |, > | 

= NGC 4286, W 
gp of 3 st, W, H 
*H 


*W,H 
nf, H 


b— Continued 


Ic 
3223 
3245 


3254 
3250 
3200 


3257 
3205 
32006 
3273 
3279 
3281 
3307 
3310 
3318 
33 IQ 
3320 
3333 
3339 
3343 
3350 
3352 
3354 
3350 
3300 
3398 
3400 
3404 
3408 
3414 
3417 
3420 
3423 
3420 
3427 
3431 
3438 
3447 
3449 
3452 
3455 
3403 
3404 
3477 
3485 
3493 
3504 
3519 
3524 
3537 
3538 
3544 
3553 


350 


il | 

nf, a defect on the plate where it 
was originally found by Frost, H; 
* * W 


NGC 4336, H 

NGC 4342 = III 95 

NGC 4341 = III 96 

NGC 4343 = Ill 94 = h1223. 
Nos. assigned according to Decl., 
in which the differences are great- 
er than in R.A., W 

nf, W, H 

*H 

NGC 4353, H 

NGC 4356, H 

1 


tou 


gi 
nf, H 
*W; = IC 3301, H 
*W,H 
*W, H 
NGC 4390, W, H 
NGC 4390, W, H 
W,H 
= NGC 4411, H 
nf, H 
k 


7 WH 

* W; found by H 

nf, W, H 

* W; found by H 

nf, W, H 

nf, H 

*H 

= HII 26?, nf NGC 4453 18/8, R 
*W; = NGC 4470, H 
7 OW, ao 

nf, W, H 

WH 

= NGC 4482, W, H 

* W; found by H 

= NGC 4402, H 

* W; found by H 

** W 

= NGC 4497, W, H 
nf, W 

“* WH 
*W 

*W 
WW 
*W, H 
nf, H 
nf, W; found by H 
or WW, 
or. A 
nf, W 

** WwW, H 
** W 














TABLE 1b—Continued 


IC IC 

3566 nf, W,H 4192 nf, W 

3569 = NGC 4561, W 4211 *W 

3572 *W,H 4338 = NGC 5334, W, Hel 
3577. ** W,H 4376 *** Hel 

3584 *W,H 4412 NGC 5594, R 
3588 * W; = NGC 4571, H 4471 NGC 5697, R 
3504 *W 4493 NGC 5747, R 


NGC 5964?, R 


3001 ~+nf, a defect on the plate where it 4551 
NGC 6014, W 


was originally found by Frost, H 4586 


bub nbu do 


3607. + * W; found by H 4013 NGC 6196, R 
3616 nf, W; = IC 3612, H 4616 NGC 6197?, R 
3647. * W; found by H 4622 nf, W 

3648 nf,H 4625 = NGC 6240, W,R 
3675 = NGC 4625, W 4626 nf, W 

3676 *W,H 4643 = NGC 6301,R 
3688 = NGC 4633, W, H 4657 nf, W 

3706 *W 4700 = NGC 6595, R, Hel 
3708 = NGC 4654, W, H 4895 = NGC 6822, W 
3712 = IC 3690, W; nf, H 5029 _ ~—snf vis, Hel 

3722 *H 5031 nf vis, Hel 

372 = I1C 3721,H 5057 *W 

3739.6 nf, H 5058 * W, Hel 

3743 *H 5061 *** W, Hel 

3700 = IC 815,H 5082 = NGC 7o10,R 
3764 = IC 816,H 5127. = NGC 7102, W,R 
3790 *W,H 5153 nf, W 

3792 +~2nf, W,H 5tss nf, W 

37907. nf, H 5stso 63* W 

3801 nf, W,H 5228 = NGC 7302, W, Hel 
3804 = NGC 4711,R 5251 *W 

3823  ** W 5260 ~—s nf, Heel 

3901 *W 5264 ~=nf, Hel 

3964 *W 5265 = IC 1459, W, Hel 
3909 «~* W, Hel 5294 nf, Hel 

4015 = NGC 4893,R 5308 ** Hel 

4016 = NGC 4893,R saux. WE 

4136 ~—nf, Hel 5313 ~—onf, Hel 

4156 nf, W 5306 ~— nf, 2" exp., L 

4190 «=o nf, W 5386 = NGC 7832,R 








NGC 48 
49 
51 

150 
178 
739 
1188 
IIgI 
I1g2 
1331 
1350 
1373 
1382 
1304 
1416T 
1422 
1455 
1489 
1518 
1521 
1040 
1744 
2053 
2077 
3143 


TABLE 2 


NEBULAE NOTED BY DREYER AS MISSING (Not INCLUDED 
IN TABLES la AND 1b) 


Dreyer IT: 


NGC 846 988 
874 1174 
Dreyer IIT: 
NGC 58 1575 
395 1592 
422 16019 
458 1689 
405 2054 
644 2203 
716 249I 
8590 2589 
866 2074 
885 2757 
1205 2860 
1327 3103 
1448 3205 
1523 3322 


Dreyer IV: 
NGC 1613 = NGC 1611 


3110 = 3122 


CORRECTIONS TO DREYER’S 
NGC 3307 


-4*, W, DIII 
ip W, D Ill 
aa; W. DI 
(2, WD an 
a9 OW, REEL 
GC 740, W 

:9, W, D III 
5.8, W, DIII 
+6/5, W, D III 
108 * p NGC 1332, W 


=F 
~§. 
5. 
7: 
2 
N 


+ 
4. 
np 
48! 
. 


3g SE 29 9 


1458 
2052 


3704 
4279 
4722 
474° 
4802 
4817 
5357 
608 2 
6200 
6294 
6328 
6308 
6403 
6450 


TABLE 3 


3h25™6, 124) 6” (1860), W 


N.P.D. 
R.A. +6: 6, W, DII 


NGC ‘ correct, WwW 
R.A. —12’, W, D III 
R.A. “ G- 
R.A. , W, DIII 
R.A. ee WW. D Ill 
R.A. —14’, W, D III 
R.A. +21’, W, DII 
R.A. —8/6, W, DIII 
N. P.D. —10’, L 
g 20™7, sad 30’, R 
N.P.D. —4‘7, W, D III 


358 


2531 
5834 


6465 
6520 
6550 
6588 
6607 
6608 
6668 
6678 
6762 
0707 
7021 
7105 
7112 
7134 


4024 
4040 
4186 
4301 
4392 
4482 
4799 
4840t 
5070 
5480 
5481 
5669 
6555. 
65798 
6580$ 
7301 
7451§ 
IC 148 

346 
1077 
1254 
2168 
5195 


N.P.D. 


3:4 ssf NGC 724 


5881 5926 
5884 7045 
7157 
7157 IC 106 
7254 130 
7520 165 
7522 453 
7055 4608 
773° 5°97 
7701 717 
7776 823 
7813 I115 
7829 1227 
1243 
1247 
1281 
1300 
1463 


POSITIONS 


N.P.D. —4‘9, W 
R.A. —7‘9, W, D 
1r® sc™, 87° 16’, 


R.A. +974, W 


ILI 
W 


12°15"3, 84° 40’, W 


7 p NGC 4380, not f,W 


RA —6.9, W, D 


Ill 


12549™2, 86° 28’, L 


NGC correct, W 
1’ sf NGC 
R.A. —21'5, W 
R.A. —23:1, W 


5077, W 


R.A. —8/1, W, D III 


1851™6, 72 
NGC correct, W 
NGC correct, W 
R.A. — 26’, W 

NGC correct, W 


N:P.D. —: 
op gett’ 
N.P.D. —30’, W 


* Value of Aa cos 6; Aa = correction to catalogue R.A. reduced to minutes of arc. 
{7 Howe’s correction (D III) to NGC not confirmed. 

t Wolf’s correction (D III) to NGC not confirmed. 

§ Bigourdan’s correction (D III) to NGC not confirmed. 


30°, W; 


—19/8, W 
ce 35™4, 108° 44’, W, D III 
(2, W, DIII 


D Il 


2, W 









































TABLE 4 


NEBULAE WHICH HAVE BEEN REPORTED Not Founp, BUT WHICH 


NGC 373 
1402 
2530 
2827 


4320 
4465 
4467 
4037 


4073 
4792 
4853 
5373 
5615 
5867 


6359 
6549 
0927 
7325 
7330 
73492 
7353 
7682 
IC 3672 


APPEAR ON Mount WILSON PLATES IN CATALOGUE POSITIONS 


** R; neb + *, W 

?a *, Hel; Sa neb, W 

Knot in s arm of NGC 2535?, R; SBc neb, W 

=IC 2460?, R; Reinmuth describes the correct neb but it is not = IC 
2460, no neb there, W 

nf by Frost, D III; SBcp neb, W 

nf by Schwassman, D III; SBa neb, W 

nf by Schwassman, D III; E4 neb, W 

= NGC 4647?, R; nf, H; the S, F neb v close f NGC 4638 was probably 
the one observed at Birr; Schwassman really measured 4638 instead of 
4637 that he records, W 

neb * 13, R; Er neb, W 

doubtful, D III; E3 neb, W 

neb *, R; neb, W 

neb *, R; Eo neb, W 

v doubtful, R; SBa neb, W 

* 14.5, not nebulous, R. Reinmuth observed a star but not the object 
found by Lord Rosse, which is 1/5 ssp of NGC 5866, W 

neb *, R; Er neb, W 

3 st, mag. 14, in line, R; SBa neb, W 

identif. doubtful, R; SBa neb, W 

? neb *, R; SBc neb, W 

nf, L; SBb neb, W 

nf, L; E2 neb, W 

nf, R; SBc neb, W 

eeF pL neb s 9’, R; SBab neb, Reinmuth refers to a group of stars, W 

= IC 809, H; So neb, W 


This work has been carried out under the general direction of Dr. 
Hubble, to whom I am indebted for advice and assistance, some- 
times in the search for an elusive photograph and sometimes in 
settling a difficult question of identity. 


CARNEGIE INSTITUTION OF WASHINGTON 
Mount WILSON OBSERVATORY 


December 1939 
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RELATIVE POPULATIONS OF 2'S AND 238 STATES 
OF HELIUM IN THE ORION NEBULA* 
ABSTRACT 


The observed equivalent width of \ 3889 of Het produced by the Orion nebula, 
together with the absence of \ 3065, leads to a lower limit for the ratio of populations 
of the 25S and 2'S states of N3/N1 => 25. This is considerably in excess of the values 
calculated by Struve and Wurm on the assumption of strict metastability for both 
levels, and thus points toward a shorter mean life for 2'S than for 25S. The matter may 
be complicated, however, by the selective absorption by the nebula of the light from 
the Trapezium stars. 

Breit and Teller’ find that, under conditions similar to those to 
be expected in a gaseous nebula, the 2'S level of He 1 should have a 
much shorter mean life than the 2S level. It is therefore of interest 
to examine briefly the observational data bearing on this point. 
The line 25S — 33P \ 3889 has been observed in absorption in all 
the stars of the Trapezium, and, from its character and displace- 
ment, it must originate in the surrounding nebula. Observations 
are most numerous for the brightest Trapezium star (#’C Orionis), 
in which the line is found to have a total absorption of 0.20 A.’ 
No lines arising from the 2'S level have been observed, although a 
careful search has been made for them. The line 2'S — 4'P \ 3965 
is well placed for observation on spectrograms taken primarily for 
d 3889, a number of which have been secured with the 32-inch 
Schmidt camera arrangement at the coudé spectrograph of the 1oo- 
inch telescope. These plates, of excellent quality, have a linear dis- 
persion of 10.2 A/mm, and it may be rather conservatively esti- 
mated that the total absorption of \ 3965 cannot exceed o.o1 A. 
Similarly, the line 2'S — 3'P A 5016 has not been found, although, 
owing to certain technical difficulties, the observations in this part 
of the spectrum are less satisfactory than those at shorter wave 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Washing- 
ton, No. 625. 

1 Ap. J., 91, 215, 1940. 

2 Pub. A.A.S., 9, 274, 1930. 
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lengths. Nevertheless, the equivalent width of \ 5016 probably can- 
not exceed 0.02 A. 

The preceding data permit the calculation of a lower limit for the 
ratio N3/N1, where N3 and Nx signify the populations of the 2S 
and 2'S states, respectively. On the linear portion of the curve of 
growth a transition having the oscillator strength f and arising from 
a level whose population is N cm~? will produce a total absorption 
A, given by 

A,« NL», 


where L is the path length through the absorbing material. (The 
geometry of the problem, as in the case of interstellar lines, justifies 
the neglect of emission into the beam.) The values of f for \ 3889 
and \ 3965 are 0.052 and 0.063, respectively. Hence a lower limit 
for N3/ Nt is given by 








N3y 0.2 | (3985) _ 0.063 
Ni ~ o.o1 3889 0.052 


This result is to be compared with the calculations of Struve and 
Wurm,‘ who derived approximate values of N3/N1 for different 
temperatures and dilution factors on the assumption that both 2'S 
and 23S were strictly metastable. Their data show that N3/N1 has 
only a rather weak dependence upon the dilution factor and does 
not differ very greatly from the Boltzmann distribution for a given 
temperature. In fact, it may be estimated that N3/N1 would not 
exceed approximately 12 and 6 for T = 104 and T = 2.5 X 104, 
respectively, even for very small dilution factors. Hence, unless the 
value of T appropriate to the Orion nebula is appreciably less than 
10,000°, the observational evidence must be considered as offering 
support to the theoretical work of Breit and Teller. The difficulty 
is that the Orion nebula is known to absorb selectively the light of 
the stars within it.’ Thus, for example, if the reddening produces an 
excess in the radiation density in the infrared, relative to that in the 


3 Hylleraas, Zs. f. Phys., 106, 395, 1937. 
4 Ap. J., 88, 84, 1938; see Table 4. 
5 Baade and Minkowski, Mt. W. Contr., No. 572; Ap. J., 86, 123, 1937. 
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blue-violet, the transition 2's—2™P \ 20582 will occur with greater 
relative frequency than is accounted for in the calculations of Struve 
and Wurm. Since 2'P connects directly with the ground state, 1'S, 
the net result will be a depopulation of 2'S, with little or no effect 
upon 2°S. Unfortunately, this possibility seems to make the astro- 
physical evidence for a shorter mean lifetime for 2'S than for 23S 


rather uncertain. 
O. C. WILSON 
CARNEGIE INSTITUTION OF WASHINGTON 
Mount WILSON OBSERVATORY 
January 1940 





NOTE ON THE POINT-CONVECTIVE MODEL 
FOR THE SUN 
ABSTRACT 

The internal temperature, density, and mass distributions of the sun have been 
calculated on the basis of the point-convective model, using tables of Strémgren and 
the Kramers’ opacity law, «=x: p/ 73-5. These tables are more extensive than Cowling’s 
and are useful for investigations into the energy production of main-sequence stars more 
massive than the sun. 

In his paper on “Energy Production in Stars,’ H. A. Bethe" quotes 
values for the central temperature and central density of the sun 
which were found by the present author on the basis of the point- 
source convective model. It seems worth while to publish this calcu- 
lation in more detail and for its own intrinsic interest. Its sig- 
nificance for a quantitative discussion of the nuclear sources of 
energy in main-sequence stars will be pointed out elsewhere.’ 

For the model under consideration the central regions must be in 
convective equilibrium, and they will therefore consist of polytropic 
cores of index m = 1.53 surrounded by “point-source envelopes.” 
In the envelope the standard equation of radiative equilibrium will 
be valid. For given M, R, L, u, and x,, the envelope equation may be 
numerically integrated inward from the surface of the star subject 


™ Phys. Rev., 55, 450, 1939. 
2H. A. Bethe and R. E. Marshak, “‘Progress Reports in Physics,” Physical Society 
of London, 1939. 
3 If we neglect radiation pressure, an assumption justifiable for the sun. 
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to the boundary conditions p = 0, T = 0, M, = M atr = R. Assum- 
ing that M = Mo, R = Ro, L = Lo, u = 1 (corresponding to the 
theoretically computed hydrogen content of about 35 per cent by 
weight), Strémgren’ has performed the integration with the follow- 
ing values of log k,: 


24.99, 24.89, 24.79, 24.69, 24.59, 24-49, 24.39. 


From his tables it is seen that for large values of k,, p>, 
T— © as r—o (the “centrally condensed”’ solutions); for the 
smaller values of x,, M, = 0 at r > o (the “collapsed” solutions). 
The physically significant temperature-density distribution cor- 
responds to some intermediate value of xk, = ky. 

To find this distribution we shall assume that the energy is being 
produced entirely within the convective core and that the outer 
radiative envelope (where the luminosity has of course been taken 
as constant and equal to the total luminosity) must be fitted to 
this convective core.° This fitting can be performed by standard 
methods.’ It is found that at the interface between the point-source 
envelope and the polytropic core we have the values 


é’ = 1.22; 03/2 = 0.776 ; RY = 1.23: 10"; 
ZT” = 15.8 X 10°; p’ = 34.7. 


6,,, is the Emden function of index 3/2 and & the “normalized” dis- 
tance. Further, at the center we have 


T. = 20.4 X 10°; Pe = 50.8. 


Also for » = 1 it is found that log, x= 24.708. The close agreement 
between these values of R’/Ro, T’/T., p’/p- and the values of the 
same quantities computed in the general case by Cowling® (by a 

4 Cf. B. Strémgren, Zs. f. Ap., 7, 22, 1933. Strémgren employed the standard model 
to compute the hydrogen content. However, it is not expected that the “point-source” 
values (at least for negligible radiation pressure) will differ appreciably from Strém- 
gren’s values; this is due to the insensitive dependence in the mass-luminosity formula 
of the average molecular weight on the energy-source concentration. 

5 These tables were made available through the kindness of Professor Gamow. 

6S. Chandrasekhar, Jntroduction to the Study of Stellar Structure, pp. 222-228, 1939. 

7 Ibid., pp. 351-355. 8 M.N., 96, 42, 1936 (cf. also n. 6). 
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somewhat different technique) is excellent. It is worthy of note that 
an application of Strémgren’s semiempirical method’ to determine 
x: leads to precisely the same result. 

Table 1 gives the temperature, density, and mass distributions for 
the sun in the region of radiative equilibrium. The values of these 
quantities in the convective core can be found simply by the use of 
tables of the Emden function of index 1.5. The superiority of Table 1 








TABLE 1 
I I ) ! r , 

t=logro (: RO O£10 logio Pp M, Mo 
o.40 12 $.077 7.862-10 I .000 
0.50 —12 6.077 8.186-10 ©.999 
0.60 —I12 6.176 8.510-10 0.998 
O90 "825 660 b es 6.276 8 .834-10 0.996 
0.80 —12 6.375 g.157-10 0.991 
0.90 —12 6.475 9.478-10 0.980 
Oe eee | 0.573 9.794-10 o.g61 
EN. iene | Ee 6.670 0.102 0.927 
0.20 —tE. 6.764 ©. 399 0.871 
O-25 =r. 6.810 0.541 0.833 
eo oe) 6.855 0.676 °.7Q91 
C86 KEE seen 6.898 0.803 0.735 
Oe SS RE ci.w 3 are 6.939 0.923 0.075 
toa) > See 6.978 1.035 ©.610 
0.50 EI 7.015 r. E37 0.540 
6.55 —11 7.051 1.228 0.470 
0.60 —II 7.083 I. 307 ©.400 
O85 7.883 1.374 0.334 
o.70 —i12 7.140 1.433 o.49 

O.75 —t1 7.165 1.483 0.221 
+ it ¢ re 7.188 I.522 0.174 
o.825—II1 7.199 1.538 0.155 




















over the Cowling table’ lies in the fact that values of 7 and p are 
given at much closer intervals in the region of radiative equilibrium 
just outside the convective core. This is important when one wishes 
to find the temperature-density distribution of more massive stars 
than the sun on the point-convective model; one may use a homology 
transformation on the sun and then correct for radiation pressure by 


means of a perturbation calculation. 
R. E. MARSHAK 


UNIVERSITY OF ROCHESTER 
9Zs.f. Ap., 2, 343, 1931; kx = x1 when @?M,/d? = o for M, = MQ/2. 
1 See Cowling, op. cit., Appendix. 
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ON THE NUMBER OF BALMER LINES 
IN EARLY-TYPE STARS 
The quantum number of the last Balmer line clearly seen on spec- 
trograms obtained with the quartz Cassegrain spectrograph of the 
McDonald Observatory,’ on Eastman Process emulsion, has been 
plotted against spectral type in Figure 1. The absolute magnitudes 
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determined by E. G. Williams? from his spectral types and luminos- 
ity groups (Mz) are shown, where available, below the name of each 
star. The absolute magnitudes of a Lyrae and ¢ Serpentis are based 


10. Struve, Ap. J., 90, 699, 1939. 7 Ap. J., 83, 333, 1936. 
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upon trigonometric parallaxes. The spectral classes are essentially 
those of the Henry Draper Catalogue, corrected in accordance with 
Struve’s sequences.’ The correlation between the number of Balmer 
lines and the absolute magnitude, within each spectral subdivision, 
is very conspicuous. Stars having appreciable rotational broadening 
were not used. 

An interesting feature of the observations is the fact that the ap- 
parent widths of the last few Balmer lines in supergiants are always 
quite small. This is conspicuous in a Cygni, but the effect is also 
present in the B2 supergiant 55 Cygni, where H.; and H., are sur- 
prisingly narrow.* On the other hand, the higher members of the 
two Het series, 2'P — n'D and 23P — n3D, are distinctly diffuse in 
55 Cygni. It would appear, therefore, that Stark broadening is ap- 
preciable in supergiant B stars. The apparent sharpness of the H 
lines must be due to the confluence of the wings, which permits us to 
discern only the central cores of fairly strong broadened lines. 

A. UNSOLD 
O. STRUVE 
McDonaLp OBSERVATORY 
AND 
YERKES OBSERVATORY 
August 1939 





NOTE ON ROSS 22 

The faint star Ross 22 was suspected to be a white dwarf (or a sub- 
dwarf) on the basis of the early spectral type, F, as determined from 
four plates‘ and from the large published motion, 0777 per year. Dr. 
Van Maanen now informs the writer that he found no evidence of 
motion for this star; the plates available at the Yerkes Observatory 
confirm this conclusion and show that the published motion was 
due to a defect. The star should therefore be dropped from the writ- 


er’s list. 
G. P. KUIPER 
March 13, 1940 
3 Ap. J., 78, 73, 1933- 
4Struve designated these lines as ‘“‘sharp”’ in his list of wave lengths of B stars 
(ibid., 90, 722, 1939). 
Ap. J., 90, 269, 1940. 





sl 














REVIEWS 


Tabelle der Hauptlinien der Liniens pektren aller Elemente. By H. KAYSER. 
2d ed. by R. Ritscut. Berlin: Julius Springer, 1939. Pp. viii+269. 
Rm. 28.50. 

Since the appearance of the first edition! of this table in 1926 many new 
measurements of wave lengths have been published, especially in the 
ultraviolet and infrared regions of the spectrum. The ultraviolet limit 
has been extended from 124 to 32.75 A, but the infrared limit of the first 
edition, 90,850 A for Nat, has been retained in the new edition. Dr. 
Ritschl points out that improvements in the photographic technique have 
produced a vast increase in the number of lines between 7000 and 
12,000 A, as well as in the ultraviolet, but that there is only a relatively 
small increase on the infrared side of 12,000 A. All wave lengths are 
given in international angstrom units. Those longer than 2000 A are given 
for air; those shorter than 2000 A are given for vacuum. Most wave 
lengths are given to o.o1 A, but there are still many for which the pre- 
cision of measurement warrants the use of only one decimal place. 

Dr. Ritschl points out that “the confusion concerning the data on 
intensities, which had already been discussed in the first edition, has 
become rather worse,”’ especially because of the tendency to increase the 
range of estimated intensities to many thousand units. ‘All intensities are 
crude estimates which have a meaning only in narrow spectral regions 
and only within a single spectrum.” Astrophysicists will find it difficult 
to understand why the experimental physicists have given so little at- 


tention to the problem of accurate intensities. 
Otto STRUVE 





Tafeln und Aufgaben sur harmonischen Analyse und Periodogramm- 
rechnung. By K.Stumprr. Berlin: Julius Springer, 1939. Pp. vii+174. 
Unbound, Rm. 36; bound, Rm. 39.60. 

This volume contains a collection of tables and numerical examples sup- 
plementing the same author’s book, Grundlagen und Methoden der 

Periodenforschung, which was reviewed in this Journal, 87, 216, 1938. 


™ See review in Ap. J., 65, 133, 1927. 
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368 REVIEWS 


Physical Meteorology. By JOHN G. ALBRIGHT. New York: Prentice-Hall, 

Inc., 1939. Pp. xxvii+392. $5.35. 

The author has presented in this book the material which he has, for a 
period of ten years, been presenting in a course in meteorology at Case 
School of Applied Science. He has emphasized the physical aspects of 
meteorology, and the book presupposes a knowledge of the principles of 
physics. One chapter on the principles of heat is included because of its 
frequent application in the study of the atmosphere. The book is divided 
into twenty chapters. The last of these, on atmospheric optics, occupies 
twenty-six pages. The book should be of interest to astronomers. 


The Concepts of the Calculus. By CARL B. Boyer. New York: Columbia 

University Press. Pp. viit+ 346. $3.75. 

The author gives a suggestive outline of the development of the basic 
concepts of calculus from their incipiency in antiquity to their final 
formulation in modern times. There are eight chapters: (i) Introduction, 
(ii) “Concepts in Antiquity,” (iii) ‘“Medieval Contributions,” (iv) ‘‘A 
Century of Anticipation,” (v) ‘““Newton and Leibniz,” (vi) ‘“The Period 
of Indecision,” (vii) ““The Rigorous Formulation,” (viii) Conclusion. 





MICROFILM SETS OF PERIODICALS 


The Committee on Scientific Aids to Learning, President Conant of 
Harvard, chairman, has made a grant to cover the cost of making a micro- 
film master-negative, on the most expensive film, of sets of volumes of 
scientific and learned journals. 

This permits the nonprofit Bibliofilm Service to supply microfilm 
copies at the sole positive-copy cost, namely, 1 cent per page for odd 
volumes, or a special rate of 4 cent per page for any properly copyable 
ten or more consecutive volumes. 

The number of pages will be estimated on request to American 
Documentation Institute, Offices of Science Service, 2101 Constitution 
Avenue, Washington, D.C. 


ERRATUM 


On page 312 of Volume go, second line of the footnote, read NaO instead of 
N20. 


